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TRANSIENT STABILITY ANALYSIS OP POWER SYSTEMS USING 
LIAPUNOV FUNCTIONS WITH IMPROVED SYNCHRONOUS 
MACHINE MODELS 


The transient stability problem arises in a 
power system when the power balance between the 
mechanical input and electrical output of CTorichronoue 
machines gets altered suddenly due to occurence of 
large disturbances such as faults, removal or addi- 
tion of large loads etc. Depending on the nature 
and location of the disturbance, the rotors of the 
machines may accelerate or decelerate. If the system 
settles down to a steady-state following such a 
disturbance or after the fault has been cleared as 
the case may be, the system is said to be stable. 
Historically this problem has been studied by the 
step-by-step methods which is still the widely used 
method with the aid of the digital computer. 

In recent years however, the Liapunov’s direct 
method has been applied to investigate the transient 



stability problem and the metbod is now acknowledged . 
a,s being able to give results in satisfactory concordance 
with those provided by the step-by-step method. Further- 
more the method has certain advantages. When the 
transient state has occured due to a fault, an impor- 
tant piece of information of interest is the critical 
fault-clearing time. The Liapunov function approach 
yields this information through a single integration of 
the faulted system differential equations while in the 
step-by-step method, repet etive integration of the 
system equations for different assumed clearing times is 
involved, each integration necessarily having to be 
carried on far beyond the fault clearing instant, even 
for the 'first swing* stability investigation. For a 
given clea.ring time the method will also indicate 
whether the system will remain stable or not. The 
liapxmov method can also be used to specify stability 
margins or stability indices, without making any actual- 
stability calculations. Because of all these reasons, 
the chances are, that Liapunov approach may prove to be 
an effective method for on-line-stability and security 
analysis. However, there are certain difficulties still 
to be overcome with the method. 1 chief dravrback of the 
Liapunov approach has been the conservative nature of 
the results. One of the contributing factors is' inherent 
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in the method itself, namely the sufficiency nature of 
the lia-punov’s theorem. Another important factor has 
heen the simplified- and therefore approximate - models 
of the synchronous machines that have "been used. It is 
■well-known that factors like damping, saliency, transfer 
conductance, flux-decay, voltage regulator and governor 
dynamics have appreciable effects on the stability of a 
power system. However, inclusion of all these factors 
in modelling the synchronous machine results in very 
complex system of equations and the, derivation of a 
Liapunov function becomes quite difficult and in some 
cases impossible as of to-day, particularly when dealing 
with multimachine systems. In the literature, many of 
these effects have not been taken into account so far 
in a comprehensive manner. This thesis seeks to cover 
some new grounds in this direction. The power system 
problem is formulated with the synchronous machine 
modelled with varying degree of details and Liapunov 
function is sought for each case. Results are obtained 
both for single manhine and nultimachine systems. 

The basic philosophy in general, is to cast the 
power system dynamics in the form of an 'indirect 
control system’ equa.tions of Lur'e's form and then 
derive a Liap'unov f'unction. Depending upon the nature 
and form of the nonlinearities, different techniques 



are called for, in arriving at a Liapunov function. 

Once the Liapunov function is .available, the stability 
region or equivalently the critical fault- cl earing time 
can be computed. The effect of parameter variation can 
also be studied with great ease and speed, using this 
liapimov function. A brief outline of the various 
chapters is given below. 

The first chapter is devoted to the statement of 
the problem of transient stability in power systems 
using the Liapunov method. The scope and objectives of 
the thesis are enumerated and organization of the thesis 
material is briefly sketched. 

The second chapter deals with the transient 
stability problem including saliency effects. A single 
machine connected to an infinite bus is considered 
first, with the synchronous machine modelled with 
damping and transient saliency taken into account. 

A Liapunov function of the 'quadratic form plus integral 
of nonlinearity’ type is developed using Kalman's 
construction procedure. This Liapunov function is used 
to determine the stability region and the critical 
clearing time. A detailed parameter analysis is 
carried out for a specific numerical example and useful 
conclusions drawn. The results with and without 
transient saliency are compared. Next a two-machine 



power system with transient saliency, is dealt with. 

Two distinct cases arise when dealing with multimachine 
systems; one is known as a system with 'uniform damping' 
while the other as a system with 'nonuniform damping' . 
Both these cases are discussed and appropriate liapunov 
functions derived. A general k-machine system with 
saliency is hriefly discussed. The difficulty in 
systematically arriving at a Liapunov function for this 
system is pointed out. It is due to the fact that when 
cast in the Lur'e form, some of the terms in the non- 
linearities cannot be determined explicitly in terms of 
the state variables. 

The third chapter considers a single machine- 
infinite bus system. The synchronous machine is at 
first modelled with the flux-decay and voltage regulator 
dynamics included in addition to damping. This model 
results in two nonlinearities requiring a matrix version 
of Ealman's construction procedure. However, a direct 
application of this procedure is difficult because of 
the fact that the nonlinearities happen to be of the 
'multi-argument' type. Using results of Desoer and V7u 
on stability of nonlinear systems with multiplicative 
nonlinearities, a systematic method is evolved to arrive 
at a Lyapunov function having the required properties 
in the region of interest. Effects of voltage regulator 


para.meters and flux decay on the critical clearing time 
are studied through a numerical example. Next a more 
general model of the synchronous machine including 
transient saliency, flux decay, voltage regulator and 
governor dynamics, is used to derive a generalized 
Liapunov function. 

Chapter four is concerned with multimachine 
system with flux decay effect included. It is observed 
that the resulting model is not amenable to the 
procedures discussed so far. Hence, recourse is taken to 
a version of 'integration by parts’ method to arrive at 
a Liapunov function. Using this technique Liapunov 
functions are derived for two-ma.chine and three-machine 
systems. The results of the two-machine and three- 
machine systems are extended to a general k-machine 
system, by induction, to arrive at an appropriate 
Liapunov function. 

Chapter five deals with power systems with 
transfer conductances. A two-machine system is 
formulated and a Liapunov function derived using 
Kalman's construction procedure for the uniformly 
damped case. 1 three-machine formulation on similar 
lines reveals that the Popov frequency criterion is 
not statisfied for this system. Indeed this is true 
for systems having more than three -machines also. 

Hence this particular line of attack has to he 
up to two-machine system. 



In chapter six a brief review of the results 
and. some concluding remarks are given. The problems 
that have surfaced during the course of this x^rork and 
which require further investiga.tion are indicated. 



CHAPTBE I 


INTRODUCTIOU 

1.1 STABILITY OB POWER SYSTEMS 

'Power system stability' is a term applied- to 

alternating current power systems, denoting a condition 

in which, the various synchronous machines in the system 

remain in synchronism or 'in step' with each other. 

Conversely 'instability' denotes a condition involving 

1 

loss of synchronism, or falling 'out of step' . The 
study of stability for purposes of analysis, has been 
divided into two major divisions; 

(i) 'Steady state stability' which is concerned with the 
stability under small perturbations or small variations 
around an operating point. 

(ii) 'Transient stability' which is concerned with the 
stability under severe disturbances such as a fault 
or sudden load changes or sudden changes in line 
reactances due to switching operations. These 
disturbances may be termed as large perturbations. 

Although stability is one phenomena, the above 
division has been made primarily as an aid to analysis - 
In fact the above two are the terminal phases of a single 
problem. For example, a test for the system for an assumed 
fault consists of two steps; A check as to whether the 
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system is ’steady state s talkie’ after tiie fault has been 
removed or isolated and the oscillations have died out, 
is first made. If it is 'steady state stable*, then it 
is tested for transient stability, generally by what is 
known as the 'first swing stability'. This consists in 
merely observing the awing curves, which are the plot of 
the rotor angles versus time, for the duration of the 
first swing and checking whether the relative rotor 
angles diverge or attain a steady state. If the rotor 
angles appear to diverge, the system is considered to be 
transiently unstable. There exists an intermediate period 
following the first swing and until the stable steady 
state condition has been reached, luring this subsequent 
period, the control devices are active and a study of the 
system behaviour during this period should consider the 
dynamics of these control devices. Such a study falls 
under the scope of 'dynamic stability'. In most cases 
the system remains stable in the subsequent swings if it is 
stable during the first swing. Therefore it is usual to 
consider a system as stable if it is 'steady state stable' 
and ’first svring stable' under transient state. A rare 
exception can however occur owing to dynamic instability 
due perhaps to abnormally slow exciter response or 

3 

inappropriate adjustment of automatic voltage regulators . 
The steady state stability investigation is straight- 
forward and relatively simple. The determination 
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tlie transieni; stability however, is more involved. In the 
most general form, it involves the combination of a solution 
of nonlinear algebraic equations describing the network and 
a solution of the nonlinear differential equations. 

1.2 IRMSIMT STiiBIlITY PROBLEM 

Transient stability can be defined as follows: 

’An electrical power system is said to be in a 
transient stability state with respect to a large 
perturbation if, after such a perturbation, it recovers 
a synchronous operating state'. (Conference Internationale 

t 

des Grands Reeoaux Blectriques a Haute Tension, 'CIGRE', 
Paris, 8-18 June, 1966). 

The large perturbation which creates the stability 
problem, may be a sudden increase in load, or a sudden 
increase in reactance of the circuit, caused for example 
by a line outage, or by a disconnection of one of two or 
more parallel lines as a normal switching operation. ■ But 
the most severe perturbation to xfhich a power system can 
be subjected to, is a short circuit, generally termed as 
’fault’. Although unsymmetricaR faults like one-line to 
ground, line to line or two-line to ground faults are more 
common in prs.ctice, a symmetrical three phase short circuit 
is the most severe of all the faults. Eor example a three 
phase short circuit on a line connecting a generator and 
a motor entirely cuts off the flow of power between the 
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machines and results in the acceleration of the generator 
and deceleration of the motor and an eventual loss of 
synchronism unless the short circuit is quickly removed. 
k similar situation arises in a multimachine system 
wherein one or more machines may experience acceleration 
or deceleration following such a fault resulting in loss 
of synchronism if the fault is not cleared. On the other 
hand, if the fault is an unsymmetrical fault involving 
only one or two lines, then some synchronizing power 
can still he transmitted past the fa,ult. In such a 
situation, under certain load conditions, the system may 
remain stable even with a sustained faiilt. However, a 
system has to he planned and designed for all possible 
contingencies including the worst conditions. Thus in 
all stability studies, a three phase short circuit must 
be considered^. 

Granting therefore, that the most severe type of 
disturbance or perturbation is a short circuit, the 
transient stability problem may be formulated for this 
eventuality in the following manner; 

Given a power system initially in steady state 
operation, assume that a fault occurs at time t^ which 
can he taken as zero without loss of generality. Two 
questions arise; 

(i) Is there a stable equilibrium position to which the 
system- settles after the fault is cleared? 
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(ii) If the answer to (i) is ’yes’, what is the maximum 

time that the fault may be allowed to remain. so that 
the system will settle to the stable equilibrium point 
after the fault i.s cleared? This is termed as the 
critical fault-clearing time or simply as critical 
clearing time. 

There are two methods to solve the above problem; 

1. The simulation method; This is presently the,, 
most widely used approach and studies the problem through 
analog simulation or by some numerical step-by-step 
integration method. Historically the network analyzer 
vms being used extensively for quite sometimes. The 
capacity of the network analyzer has been limited however, 
thereby restricting its use for reasonably small systems 
having limited number of generators. But with the growth 
of large power systems due to addition of generators and 
increasing interconnections of various existing systems, 
power system engineers have changed over to the use of 
digital computers for stability analysis by munerical 
method. The standard numerical integration method 
approach divides the transient stability problem into 
two phases; 

Phase A: This phase consists the study of the evolution 
of the system from the instant of fault to an 
assumed clearing time t^. 
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Phase B; This involves the study of the evolution of the 
S 3 rstem in its post-faiilt period. 

The model can have detailed generator representation 
to include flux decay, voltage regulator, governor action 
etc. For an assumed clearing time t^, the transient 
stahility is examined. Clearly, if t is less than the 
critical clearing time t , the system will remain stable.. 
In general however, the critical cl ea, ring time t will 
not be kno-vm apriori. It is determined by several trials 
of assumed values of clearing time Thus the estima- 

tion of t^ involves repetitive computations because of 
phase B and hence time consuming. 

2, ’Direct analysis' methods By this method, the 

stability is investigated without actually solving the 

1 2 

differential equations. The equal area criterion ^ , 
the phase-plane method the energy integral method^ 
and now the Liapunov's direct method " fall under this 
category. While the first two of the above are applicable 
only for second order systems, the last two methods can 
be used for higher order systems and thus for multimachine 
systems. 

1.5 LIAPUNOV’S DIRECT METHOD 

The Liapunov method is more general than the 
energy integral method. In fact, energy integral 
happens to be a particular case of a specific Liapunov 
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function, viz. >?'hen the total system energy is chosen as 
the Liapunov function. But in general, Liapunov function 
will not he restricted to the energy function. - Indeed,: 
the more acciprate Liapunov fianctions which yield larger 
stability regions, are not ezactly the energy integrals 
of the system. It can therefore he safely speculated 
that among all direct methods, Liapunov approach is the 
only method, that can possibly compete with the standard 
step-by-step method. There are good reasons to visualize 
such a prospect as discussed below. 

The. Liapunov’s direct method eliminates the phase B 
computation discussed in method 1 above. Instead, - a 
st8,bility criterion is used.. This criterion is based on 
the construction of a Liapunov function V for phase B 
system around the post-fault equilibrium point and its 
comparison with a limit value 0. 

let the dynamic equations for the faulted and 
post-faxxlt condition be denoted respectively by 


z = B^(x) 

(A) 

i = I'2(x) 

(B) 


At the outset, a Liapunov function Y{x) is constructed 
for system (3) and a constant C is computed such that 
V(x) < C constitutes a region of attraction around the 
post-fault stable equilibrium point. The constant C is 
generally found by evaluating the function V(z) for 
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system (B) at the unstable eq.uilibrium point nearest to 
the stable post-fault equilibrium point. Fow, given a 
clea^ring time t^, the faulted system (A) is numerically 
integrated and the value of the state vector x is 
computed at the end of t . Fow with this value of the 
state vector, the liapunov function V is evaluated. 

This value of V is compared with C. The system is 
considered stable if V <. G. The critical clearing time 
also can be estimated by evaluating the V-function at 
the end of each integration step and comparing it with 
C. That value of t at which 'V(x) = 0 gives the 
critical clearing time t^. The computation of t^ thus 
involves just one integration of system (a) and does 
not require any integration of system (B) unlike in the 
method 1. The saving in computation time 'can be 
significant, particularly if a large system is to be 
examined for its transient stability property for various 
contingencies. Here the role of the Liapunov method is 
to be considered as complementing rather than supplementing 
the actual detailed simulation. A preliminary screening 
by Liapunov method may reduce the number of studies to be 
done in detail. Further, the influence of the various 
parameters can be thoroughly investigated with great ease 
and speed. Another advantage of the Liapunov method is 
that it can be used to specify stability margins or 
stability indices without making actual-stability 
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calculations. Gradually it is being acknowledged that 
the Liapunov method is able to give results in satis- 
factory concordance with those provided by the step-by- 
step method. The indications at present are, that this 
approach may still succeed in being accepted as an 
effective tool for on-line-stability and security 
analysis. However, there are certain difficulties 
still to be overcome with this method. 

A chief drawback of the Liapunov method has been 
the conservative nature of the results. One of the 
contributing factors is inherent in the method itself, 
namely the sufficiency nature of the Liapunov theorem. 
Another important factor has been the simplified- and 
therefore approximate-models of the. synchronous machines 
that have been used. It is well-known that factors like 
damping, transient saliency, transfer conductance, flux, 
decay, voltage regulator and governor dynamics, have 
appreciable effects on the stability of a power system* 
But, inclusion of all these factors in modelling the 
synchronous machine results in a set. of system eq. nations 
which are highly complex and the derivation of a 
reasonably good Liapunov function becomes quite difficult 
and in some cases impossible as of to-day. This is more 
so in multimachine system. Many of these effects have 
not been taken into account in a comprehensive manner 
in the literature while applying the Liapunov's method. 
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This thesis seeks to cover some new grounds in this 
direction'. The power system problem is formulated with 
the synchronous machine modelled with varying degree of 
details and a Liapunov function sought for each case'. 

Both single and multimachine systems are considered, 

Before giving an outline of thesis, some preliminaries 
are discussed. 

1.4 THE COBSTBUCTION OP LIAPUNOV FUNCTION 

Starting with the pioneering work of El-Ahiad and 
7 S 

Nagappan and that of Gless , there have been a number of 

researchers vrho have developed or proposed various Liapunov 

functions for studying transient stability of power systems, 

although many of the Liapunov functions either reduce to or 

are simplified versions of the one developed by Bl-Abiad 

7 

and Nagappan . Various techniques have been used by 

different authors to arrive at the Liapunov functions. 

Trial and error and energy consideration have been used 

in references 7, 8, 9, 12 and 25* References 10 and 15 

29 18, I’J ■ 

use variable gradient method . Dharma Rao ^ has used 

Cartwright method, Aizerman method, method of integration 

by parts^^ and method of Infante and Clark^^. Zubov's 

method-^ has been used in references 11 and 26, Recent 

58 

interest has been to use Popov’s frequency criterion 
as a necessa.ry and sufficient condition for the existence 
of a special type of Liapunov function of a quadratic form 
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34 35 

plus an integral of the nonlinearities ’ . Based on 

34 

this criterion, Kalman's construction procedure"^ has 

been successfully used to construct Liapunov functions 

for single machine and multimachine systems by Pai 
20—22 

et.al . In a parallel development, the results of 

33 

Moore and Anderson^ for the stability of multinonlinear 
systems have been used to get Liapunov functions by 
¥iliems^‘^~^'^ and Pai and Murthy^^’^^ for multimachine 
systems. The methods that yield better Liapunov 
functions for single machine system suffer from the 
limitations that it is very difficult to extend these 
methods to multimachine systems or systems with more . 
detailed representation of synchronous machines. As 
this thesis will demonstrate, when more detailed machine 
models are used, it becomes necessary to modify some of 
the known methods, to suit the particular problem. 

1.5 SOME GEHERALITIES 

Consider a nth order dynamic process which can 
be represented in an n-dimensional state space as: 

X = A X + b f (o) 

CT = c^ X (1.1) 

where A is a nxn real constant ma.trix, b and c real 

constant n-vectors and x is the n-state vector, c is 
a scalar. f(o), the nonlinearity which is also called 
'characteristic' , is an arbitrary, single-valued 
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piecewise continuous real function, defined for all 
real values of a and satisfying 

f(e) = 0 

0 < k ■ (1.2) 

where k is either a finite positive number or infinity. 

In the latter case it is equiva.lent to the inequality 
of (o') ho. It ■will be further assumed that the 
system (l.l) satisfies the usual conditions necessary 
to ensure the existence and uniqueness of a solution 
for all t h 0, starting from a.ny initial state x{0) = x°. 

If in (l.l), the matrix A is strictly Hur'witz, 
by vjhich is meant that the roots of the characteristic 
equation 

Isl - A| =0 (1.3) 

all lie on to the left of the imaginary axis of the 
complex s-plane, the system is kno-wn as *the principal 
case'^^ or sometimes called 'direct control' system^"^. 
Suppose the characteristic equation (1.3) has one zero 
root. It is then possible by means of a linear 
nonsingular tra.ns format! on to transform the original 
system (l.l) into a system 

y = F 2; - s f(c^) 

5 = - f(a) 
or = h^ y + P ? 


(1,4) 
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where now P is a (n-l)x(n-l) real constant matrix which 
is Hurwitz, i.e. 

! si - pI = 0 (1.5) 

has all its roots on the left half of s~plane. g and h 
are real (n-1 )-vectors, y is a (n-l)-state vector and 
C the nth state variable. When system (l.l) is such 
that the characteristic equation (1.3) bas one zero 
root and tho remaining (n-l) roots on the left half of 

s-plane, the system is known as the 'simplest particular 

'36 37 

case'-^ or 'indirect control' . In general, if the 

system (1. 1) is such that some of the characteristic 

roots of (1.3) are on the imaginary axis (including 

zero roots) and the rest in the left half of the s-plane, 

it is called a particular case of (l.l) or briefly 

'particiilar case'. Por particular cases, the class of 

admissible nonlinear functions have to be narrowed 

down to the inequality 

0 < — h (a ^ 0) if k is finite, 

and o'f(cr) ^ 0 (c ^ 0) if k is infinite 

( 1 . 6 ) 

The systems (l.l) and (1.4) are equivalent as long as 
the transformation is linear and nonsingular and therefore 
the stability properties of (I.4) are identical to those 
of (1,1)^^"^^. 
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l.,6 POV/BR SYSTEM PROBLEM 

A single machine infinite bus system under 
transient condition can first be formulated in the form 
of (l.l) for phase B referred in section 1.-2.- It will 
be observed, that the power system problem belongs to 
the particular case of (l.l). By suitable transformation 
the problem will be represented in the form of (1.4)* 

In general, a power system formulation may contain more 
than one nonlinearity either because of more detailed 
representation of synchronous machine or because of the 
presence of more number of machines. The most general 
form of the system equations would then be 

X = A X + B f ( 2 ) 

0 = X (1.7) 

where 0 and f(o;) are now m-vectors, and B and C are 
nxm matrices. Furthermore, if the power system has k 
synchronous machines (called hereafter a k-machine 
system), one would discover, that the matrix A in (1.7) 
is such that the characteristic equation (1.3) will now 
have k zero roots and the rest of the roots on the left 
half of s-plane. It is still possible to transform this 
most general system into what may now be called the 
matrix version of (1.4) as 

y = F y - G £( 0 ) 

1 = -N £( 0 ) 

T r 

= H-" y + p i 


a 


(1.8) 
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•where F is now a (n-k)x(n-k) real constant strictly 
Hiirwitz matrix, G, N, H and P are r eal > cons tan t » 

(n-'k)xm, kxm, (n-k)xm and mxk matrices respecti'wely. 
y and £ are (n-k) and k-irectors respectively and f 
and 0 being m-vectors. For convenience we will call 
( 1 . 4 ) and ( 1 . 8 ) as the 'standard forms'.. 

The differential equations for the phase B are 
cast in the form of (1.4) or (1.8) as the case may be-. 

We now seek a Liapunov function consisting of a 'quadratic 
form plus an integral of the nonlinearities* type. The 
relevant Liapunov theorem on which we base our stability 
analysis of power system is given in Appendix A. It is 
appropria^te to note that the nonlinearity? .f(o) in a' 
po'!»?er system problem, does not satisfy the 'sector 
condition' (1.6) in the whole state space. However, 
around the post-fault equilibrium point which will be 
taken as the origin of the state space, there is a 
region in which the condition (1.6) is satisfied and 
therefore the stability properties can be studied by 
the Liapunov method in this region. Indeed, this region 
of stability will be of interest since the critical 
clearing time is computed with the aid of this. 

1.7 SCOPE OF THE THESIS 

A brief account of the work that is reported in 
the following pages is outlined below. 
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The synchronous machine model including transient 
saliency is dealt with in the second chapter. A single 
machine connected to an infinite bus is considered first, 

A Liapunov function is developed using Kalman's construc- 
tion procedure^^. This Liapunov function is used to 
determine the stability region and the critical clearing 
time. A detailed parameter analysis is carried out for a 
numerica,! example. For the sake of comparison a similar 
parameter analysis is also made with transient saliency 
neglected, A two-machine system with both machines 
modelled with transient saliency considered, is studied 
next. When dealing with a multimachine system, two 
distinct cases arise; one in which the damping is 
'uniform' , by which is meant the damping ratio 
is the same for all machines; and the other case wherein 
the damping is ' nonuniform ' , meaning thereby that the 
damping ratios for different machines are different.. 

The two-ma,chine system referred to above, is studied for 
both cases. A general k-machine system is then briefly 
discussed., indicating the difficulties involved in 
systematica, lly arriving at a Liapunov function for this 
system. 

A single machine system with greater details in 
machine representation, is the theme of the third chapter,. 
First a model with flux-decay effect and voltage reulator 
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action included, is investigated This prohlem turns out 
to be mult in onl in ear system. Further the nonlinearities 
happen to be of the 'multi-argument' type.- Therefore, a 
direct application of a matrix version of Kalman's 

•tQ'zq 

construction, procedure^ ’ is not possible. However, a 
modification of this procedure with the incorporation 
of results from Desoer and Wu"^*^ on stability of nonlinear 
systems, provides a systematic method to arrive at a 
Liapunov function. With the Liapunov function developed, 
the effects of voltage regulator parameters and flux— decay 
on the critical cleering time are illustrated through a 
numerical exa.mple. Next a more complete model of a 
synchronous machine to include transient saliency,- 
flux-decay, voltage regulator and governor dynamics,- 
is used to derive a generalized Liapunov function. This 
Liapunov function naturally subsumes the previously 
derived Liapunov functions as special cases. 

Chapter four is concerned with the multimachine 
system with flux-decay effect included. The resulting 
model is not amenable to the procedures discussed in the 
previous chapters. A version of 'integration by parts' 
method is found suitable. Using this technique, Liapunov 
functions are derived for two-machine and three-machine 
systems. These results are extended by induction to a 
general k-machine system. 
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Chapter five deals with a power system with 
transfer conductance, A two-machine system is formulated 
and a Liapunov function derived through the matrix 
version of Kalman's construction procedure, A three- 
machine system is then formiiLated. 

A brief review of the results, .identification 
fjf some unsolved problems and suggestions for future 
work form the sixth and the concluding chapter. 



CHAPTER II 


BPPBCT OP SALIEHCY OH TRANSIEMT STABILITY 

2.1 IITRODUCTIOH 

In most studies of transient stability, a simple 
model of the synchronous machine represented by constant 
voltage back of a transient reactance is used. This 
reactance is assumed constant, regardless of the angular 
position of the armature reaction with respect to the 
rotor. This is equivalent to assuming that the direct 
axis and quadrature axis transient reactances are equal 
(x^ ~ that a slowly decaying component of rotor 

flux linkages exists in the direct as well as the 
quadrature axis. In reality however, most of the rotor 
circuit flux linkages decay relatively rapidly except 
those supported by the main field winding which has a 
relatively large time constant. It is therefore more 
correct that only the main field winding flux linkages 
are held constant during the transient period and the 
other more rapidly decaying flux linkages in the rotor are 
neglected. The synchronous machine may then be said to 
have' transient saliency’ (x^ ^ x^)^. It is to be noted 
that transient saliency is present in both salient pole 
and cylindrical rotor machines, although the latter has 
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practically no synchronous saliency. The synchronous 
saliency is reflected in the electrical output in the 
form of an additional term for the salient pole machine. 
¥e will consider salient pole machines in the subsequent 
sections. 

The effect of saliency has been discussed in 

greet detail in earlier works, notably in references 2 

and 41. Many others have used salient-pole model in 

their study of transient stability^^''^^’^^^^'^. The aim 

of the present chapter is to develop a systematic 

procedure for constructing a Liapunov function for a 

system model including transient saliency. Liapunov 

functions are developed for a single machine-infinite 

bus system and a two-machine system using Kalman's 

34 

construction procedure . Por the latter both uniform 

and nonuniform damping are considered. These Liapunov 

43 

functions can be used to find regions of stability . 

The critical fault clearing time for each system is 
computed by numerical integration of the faulted system 
equations until the trajectory reaches the boundary of 
the respective stability region. A detailed parameter 
analysis is carried out for the single machine case. 

The complexity of the problem when extended to a 
k-machine system with transient saliency is indicated. 
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2.2 SIIGLB MACHINE COHKBCTED TO IIEINITE BUS 
2.2.1 P omul at ion of the Problems 

The dynamic eq.uation the single machine system 
can he written as 


M6+I)6 = P, 

I e 


( 2 . 1 ) 


where P^, the mechanical input>is assumed to remain 
constant. The damping- cc<='fficiont' ia'als ©'•assumed ‘tfl 
be.': constant ; Variable damping can be taken into account 
as discussed in reference 41 at the expense of more 


analytical computations 


26,27 


. The electrical po'wer output 


P of a salient pole machine connected to an infinite bus 


can be computed from' 


44,2 


E e. 

P = sin a... + — sin(6-a-,-) 

e 11 


( 2 . 2 ) 


where Eq is the voltage behind the quadrature axis 
synchronous reactance and e^ is the infinite bus voltage. 
Zp! and Zp2 driving point and transfer impedance 

respectively, with quadrature axis synchronous reactance; 


(Y — E _ Q 

“ll 2 ^11 ^ 


a. . = 


12 " "'ll 


©, T and ©-, o are the 


impedance angles of Zpp and Zp^ respectively (radians). 

The voltage Eh is related to a voltage E* which is the 
y q 

quadrature axis component of the voltage behind the 
direct axis transient reactance (E’ is a voltage 

q. 

considered proportional to field flux linkages) by the 
2 

relation 



22 


(z'-x 

E' = [1 + — V- 


. = [1 ^ 


•. “'“.p 

(2.3) 


If resistances are neglected, the electrical power in 

2 

terms of E' and e, can be shown to be (Refer to the 
q t 

phasor diagram of Figure 2.1). 


E' © 4 . e7(xi-z ) 

_l_i sin 6 + .i -i -S sin 26 

2(x^4-2g)(x^+Xg) 


= sin 6 + F 2 


wrhere P-, = 


E' e 




; Pj = 


2(x^+X^)(x’+Xg'> 


(2.4) 


Here z is the external rn.actance through which the 

0 

machine is connected to the infinite bus. It should be 
noted that P in (2.1) has to be evaluated with P, and 
P^ computed for the particular condition of operation 
such as prefault, post-fault or faulted condition, with 

the corresponding network impedances. The equilibrium 

** ♦ 

points of the post-fault system where 6=6=0, can be 
determined by solving for 6 from 

Pj - P^ sin 6 - P 2 sin 26 = 0 (2.5) 

s 

let 6 be the value of 6 at the p*st-fault stable 


equilibrium point. 
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Defining 

X = 6 - 6®^ (a, 6) 

the equation (2,1) can he rewritten as 

,M X + D X = Pj ” sin(x+6®) - P^ sin(2x+26®) 

(2.7) 

Fow, define the state variables x^ and x^ as 

x^ = X 

x^ “*■ X ,(2.8) 

let cr = 

P P P 

f(a) = ^ sin(a+6®) + ^ sin 2(a+6®) - ,..(2.9) 


Using (2.8) and (2.9), the system equations in the. state 
space can he put as 


.. • 1 

— — ; 




P” "T 

i 

0 1 


^1 


0 

j = 




4” 



0 -D /M 


^2 


-1 


c 


- [1 



( 2 . 10 ) 


This system can he represented hy the block diagram shown 
in Figure 2.2. The nonlinear function defined hy. (2.9) 
satisfies the sector condition f(0) = 0 and af(o) > 0 
for a ^ 0 for the range (see Figure 2,3a and 2.31) 
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q-axis 


FIG-. 2.1 PHASOF DIAGRAM OF A SAIIMTT POLE' 
GENERATOR. 



FIG. 2.2 SYSTEM >iITH A SINGLE NONLINEARITY. 
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^(2ix-g’^-6®) < g < (6^-6®) (2.11) 

Popov’s criterion^^ can therefore he applied for analysing 
the system in the region where the sector condition is 
satisfied hy the nonlinearity. Indeed the operating region 
of our interest falls within this region. System (2,»10) 
is similar in form to (l.l). We transform this into the 
form (1.4) hjr defining a new state variable 


D 

-- X. + 

Ml 2 


( 2 . 12 ) 


V/ith the introduction of this state variable and 
eliminating z^ from (2.10) , the system equations reduce 
to 


^2 - 


£ X 

M 2 


f(c?) 


? = — f ( 0 ) 

MMp /'Ol'ai'i 

0 = - ^ x^ + \ 2 , 13 } 


A Liapunov function is now to he derived for the system 
( 2 . 13 ) which is in the ’indirect control' form of the 
lur’e sj^stem , 


2.2.2 Liapunov Function for the Lur'e Problem with a 
Single honlinearity: 

rz ft 

Consider the Lur'e problem defined by^ 

X = F X - g f ( 0 ) 

I = -f(0) 

T 

0 = h X + P ? 


(2,14) 
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where F is a real nxn strictly Hurwitz matriz, x, £ and 
h are real n-vectors and a,? and P are real scalars. 
f( 0 ) is a real valued continuous function which satisfies 
a sector condition 

f(0) 0 and o f(0) >0 a ^ 0 (2.15) 

For the system (2.14) Popov^^ enunciated the following 
theorems 

'Assume F is stable and that P > 0. Then (2.14) 
is globally asymptotically stable if the condition 

Re { ( 2a D + j “ P ) [ h^ ( j I-F )'"^g + 

— J CO . 

for all real “ (2.16) 

holds for 2a P= 1 and some p R-Ol* Differentiating the 
third equation in (2.14), a system equivalent to (2.14) 
tan be written as 

X = F z - £ f(p) 

a = h^ F z - (h^ g +p ■) f(0) (2.14*^ 

Popov showed that for the system (2.14*) (or equivalently 
for (2.14)) the most general form of a Liapunov function 
which is of the kind described as 'quadratic form plus 
the integral of the nonlinearity', takes the form 

Y(z, o) = z^ L z + a(0-h^ z)^ + p / f (o) da (2.17) 

0 
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He further established that the freq.uency condition (2.16) 
is necessary for the existence of a liapunov function of 

. , 34 

the type given in (2.17). Kalman . subsequently showed that 
(2.16) Is also sufficient for the existence of the above 
type of V-function. Furthermore, utilizing and extending 
the results of Popov^^, Yakubovitch^^ and laSalle^^, 

Kalman developed a systematic procedure for constructing 
the above type of Liapunov function. The procedure which 

relies on what has come to. be known as Kalman-Yakubovitch 

47 

lemma can be summarized as follows : 

Define Z(s) = (2a P4- Ps) ¥(s) (2.18) 

where W(s), which is the transfer function of the linear 
part of the system relating o and f(a) (see Figure 2.2) 
is given by 

¥(s) = h^ (si - F)“^ g + ~ P (2.19) 

The frequency condition (2.16) is now equivalent to the 
f oil owing s 

) + Z(-o w)] h 0 for all realm (2.20) 

If now for the system (2.14), the above frequency 
condition is satisfied, then a Liapunov function exists. 

This Liapunov function can be constructed in the following 
steps: 

(i) Factorize the sum -^[ 2 ( 3 ) + Z(— s)]in the form 

•|■[Z(s) + Z(-s)] = T(-s) T(s) (2.21) 
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(ii) Determine the scalar y from 
Y = P(p + h^ g) 


( 2 . 22 ) 


(iii) Solve for the elements of n-vector u from the 
identity 

rp —1 

l(s) - Yy = (si - P) g (2.23) 

(iv) Now solve the liapunoxj" matrix equation 

F^L + 1 P = - u (2.24) 

for the syminetric positive semidefinite (possibly positive 
definite) nxn matrix L. Knowing 1, the Liapunov function 
is now obtained as in (2.17). The Liapunov function thus 
obtained will satisfy the sign definite properties (V ^ 0 
and Y "^0) required for the absolute stability provided 
the nonlinearity satisfies the sector condition (2.15). 
Zalman imposes the requirement that (P, g) be completely 
controllable and (P, h ) be completely observable. These 
restrictions have been relaxed subsequently by Meyer^®.’ 


2 , 2.3 Application to the Problem in Section 2.2,1: 

The power system problem described by (2.13) is 
seen to be identical to the lur*e system (2,14), with 
the following identities: 


P = 


D . 

M ' 





(2.25) 
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¥ith these identities, the transfer fimction W(s) given 
hy (2.19) for the linear part is seen to he 


¥(s) 


1 

s(s + A) 


vrhere 



Choosing 


a =..D/2M 


( 2 . 26 ) 


(2,. 27) 


the left hand side of the freq^uency condition (2.20) 
becomes 

i[Z(j )] = (2.28) 

Thus the frequency condition (2.20) will he satisfied 
if >^811 (2.29) 


p can he so chosen as to satisfy this condition. 

how the Liapunov function 7 is constructed: 

(i) ¥ith Z(s) as in (2.19) > it can he verified that 

t[2(s) + Z(— s)] = ^ 2 ' 2\' 

( ~ S 4- ^ ) 

Factorizing as in (2.21), we have 

T(s) = f(^ P - l)/(s + ^ ) (2.30) 

(ii) Using (2.22) it can he seen 


Y = 0 


(2.31) 
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(ili) Noting that in this problem u is a scalar^ we have 
on solving (2,33) 

u = - V'( - 1) (2.52) 

(iv) Again 1 is a scalar now (say ). The liapunov 
matrix equation (2.24) now becomes a scalar 
equation and solving this we have 


i±L^i l 
2 A 


(2.53) 


Hence the liapunov function as given by (2.17) is 


/Ad-i\ . , Mx~a <y 

(j, a) = z" 4- ^(o + -^) + p / f(a)da (2.54) 


V( 


■ 2 ''^' “ 2 ^ "2 ' 2 ^" " 1 


This is the most general form of liapunov function with 
P — 1/ ^ . For the particular choice of p = 1/ ^ we have 


-p. M X^ 2 M ^ 

Y(z^,a) = g f(a)dcT (2.35) 

¥ith f(a) as in (2.9) » (2.35) becomes 

M X 2 ■ P 

V(x 2 ,a) = + |C||^(cos6®-cos(a+6®)) + 

P P 

2 |p(cos 26®-cos 2(a+6®)) ~ o’] (2.36) 


For the power system model under consideration, the 
nonlinearity violates the sector condition as it leaves 
the first and the third quadrants at values of o equal 
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to cr, and a, respectively (Figure 2,31) • such a 

ix 

system there is only a finite region of asymptotic 
stability around the origin. This region R around the 
origin is defined by the inequality 

V < c ■ (2.37) 

where C = min (C^^, C2) 

with Cn= min V(x2> or) 

0=0^ 

C2 = min Y(x^, a) (2.38) 

43 

Determination of and using lagrange multipliers 

20 

has been discussed by Pai et.al. . For the system under 

47 

consideration, ^ and it can be shown that 

O' £ 

C = I / f (o)da (2.39) 

where 0£ = 6'*'^ - 6^. 6'^ is the value of 6 at the post- 
fault unstable equilibrium point. Using (2.9) and (2.39) » 
we have 

P P P 

Cj = cos6®-cos6^] + ^(cos 26®-cos 26'^) - ^(<5'^-6®) 

(2.40) 

The region of asymptotic stability is given by 
^(^o+z^)^ + ^[P^(cos6®-cos(o+6®) ) + 

o o 

■^(cos26®-cos2 (0+6° ) ) - Pj o] < 


( 2 . 41 ) 
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Por a more general clioice of p satisfying p ^ 1/ ^ , 

corresponding Liapunov function can be obtained. 

Critical Clearing Time: 

The prefault operating point 6^ is first computed 
for the initial operating condition from load flow,. The 
stable and unstable post-fault equilibrium points are 

7 

determined by the well-established methods such as 
steepest descent method, fletcher Powell or Pletcher 
Reeves niethod. Then, starting from this prefault 
equilibrium state, the system differential equation (2.7) 
is numerically integrated for the faulted condition till 
the boundary of the inequality (2.41) is reached. This 
gives the critical clearing time* 

Parameter Variation: 

Equation (2.7) corresponds to a particular set of 
parameters. While the system parameters M and D are 
decided upon when designing, and nre dependent on 
the initial excitation and location of the fault while 
Pj depends on the load demand, at the instant of fault. 

In an overall system planning and study, it may be of 
interest to investigate how variations in these 
parameters reflect on the critical clearing time. Such 
a parameter analysis is illustrated through a numerical 
example in the next section- The influence of other 
par8.meters such as the various reactances, which also 
affect stability, can also be similarly investigated. 
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liumerical Example: 

Figure 2.4 shows a power system with a three phase 
fault on one of the two transmission lines . Each 
transmission line has a reactance of 0.2 p.u. The 25 MYA 
generator has the following constants: 

H =; 2.76 MJ/MVA ; x^ = 1.0 ; x^ = 0.65 ; = 0.5 ; 

Frequency = 60 Hz. 

The transient stability of the above system is investigated 

for various situations. On the lines indicated in the 

49 

previous section, a parameter analysis is made. For an 

assumed three phase to ground fault at the middle of one 

of the transmission lines, the values of D, M and Pj- are 

varied. Similarly keeping a set of fixed values for D, 

M and Pj, the fault location is varied. For each set of. 

parameters Y{x^f or) and are computed. The effects of 

all these variations on the critical clearing time t are 

c 

investigated as shown in Figures 2.5 to 2.8. In studying 
the effect of fault location, a parameter a is introduced 
for convenience, a is a number between 0 and 1 and is 
defined as the ratio of the amplitude of the fundamental 
component of the electrical power during the faulted 
state to that during the prefault condition. Thus a in 
some way, is a measure of the reduction in the power output 
due to the fault which in turn, is dependent on the 
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fault 



FIG. 2 .^ A SYNCHRONOUS MACHINE -■ 
INFINITE BUS SYSTEM. 





FIG. 2*6 EFFECT OF INERTIA ON 



LEGEND 

Saliency 
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location of the faTolt. Por comparison, the critical 
clearing time for the above conditions, with transient 
saliehcy neglected, are also computed and plotted. 

Discussion on the Results: 

The parameter analysis conducted as described 

above, yields very useful information. Figures ,2.5 to 

2.8 show clearly, the effect of the variation of the 

various parameters on the critical clearing time t^. 

Saliency in general has an effect of increasing critical 

clearing time. The extent of increase however, depends 

on the other parameters also. The effect of damping far 

instance is also to increase t . But this effect is more 

c 

pronounced as. the input power Pj is reduced or equivalently 
if the initial load on the machine is less. The results 
indicate that if the machine operates with outputs less 
than 1 p.u. , it may be worthwhile to consider introducing . 
higher damping so as to increase the critical clearing 
time. However, if the machine is required to operate at 
1 p.u. or higher output power, an increase in damping 
coefficient D does not seem to significantly affect t^ 
(Figure 2.5)* These observations may be difficult to 
translate in terms of physical parameters since damping 
is known to be not constant in general. Nevertheless it 
offers some insight into the otherwise highly complicated 
phenomena. Figure 2.6 shows how t varies with inertia 
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constant H, ■whicii is linearly related to M, appearing in 
tile differential equation. A higher value of inertia 
constant results in a higher value of t^. This is 
consistent with physical considerations; for, a machine 
with higher inertia constant will take more time to reach 
the critical clearing angle, implying thereby a greater t . 
Figure 2.7 shows effect of P_ on t with H, .I and a held ' 

-L C 

constant. As can be expected, lower the value of P^, the 

greater is the vs^lue of t . Here it is observed that the 

o 

effect of saliency is more pronounced at lower value of 
Pj. Thus, at lower loads on the generator, saliency 
seems to increase t considerably. The effect of varying 

Co 

the location of fault - which corresponds to varying a - 
with H, D and P^ held constant, is shown in Figure 2.8. 

It is observed that higher value of a (i.e. less 
reduction in electric power output during the faulted 
condition or which is equivalent to a fault farther away 
from the buses), greater is the value of t^. This is 
understandable, because, the farther the fault away from 
the buses, less severe it is. It can also be seen from 
the graphs, that saliency has greater influence on t at 
higher values of a. In practice Pj is dictated by the 
prefaixLt load demand and a signifies the fault location. 
These parameters cannot be pr eassigued . However, one 
can prepare a chart for ready reference predicting the 
behaviour of the system for various contingencies. On 
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the other hand, inertia constant and damping and to some 

degree the relevant reactances are controllahle parameters. 

It may therefore te worthwhile to make a study like the 

ahove during planning and designing a system.. Some of 

2 

the above conclusions have been reported in Crary 
through extensive simulation. However, Liapunov method 
in conjunction with digital computer affords a quick way 
of obtaining extensive results on effects of parameter 
variation. 

2,3 TWO-MACHIHE SYSTEM 
2.3»1 Formulation; 

The theory discussed in the foregoing sections 
can be extended to the two-machine system, with both 
machines having saliency-. The equations of motion for 
the two machines can be written as 


M. 

1 


6, + E. 6^ = - P 


mi 


ei 


(i=l,2) 


(2.42) 


where the subscript i stands for the ith machine. Using 

these subscripts for voltages, currents and reactances of 

2 

the corresponding machines, it can be shown .that the - 
power output of machines 1 and 2 will be given respectively 
by 
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= ^7 sin “'ll + ^ sirL(6^-62'-a3_2) 


el 11 Z^2 


e2 


^ sin (X22 + sin(62-6^-a^2> 


(2.43) 


Eqi and s-ne the voltages hack of the quadrature axis 

synchronous reactances of machines 1 and 2 respectively. 

(i=l,2) and ^ne the driving point and transfer 
impedance respectively. It is to be noted here that 
Z . . = Zj. and also a. . = “tt* and 6^ are the angles 

Ij J. 2 . Ij jX X xc 

of the voltages and respectively from the 
synchronously rotating reference axis. The voltages 
E^^ and are to he determined from the following 

two equations: 


E^l = a^^ Eq^ '^^^ 2.2 sin (6^2 + ®12^^®Q2 


® q 2 ^"21 ® Q 1 ■*' ® 12 ^ ^-22 ® Q 2 


(2.44) 


where E' and E’ pj the voltages proportional to field 
flux linkages of the machines 1 and 2 respectively, 
are known and held constant (hy the constant flux linkage 
theorem). Eurther 6. . =* 6. - 6 .. The coefficients are 

X J X j 

given hy 


ail 


^^dl - 


Z 


11 


sin j a^2 


(x -xIt ) 
^ ql dl 

^12 


8, 


21 


= ~ ^d2) 

Z 12 


(x’ -X^,) 

a^^ = 1+ - 77- " — sin Q 


'22 


'22 


22 


(2.45) 



44 


Neglecting -the resistances for ihe sake of simplicity, 
we have = 4>^2 “ '®12 “ ^11 “ ^22 ~ °^12 ~ 

Then equation (2.43) reduces to 


E--, E^o 

P =: g 5 ^ P 

el Z ^2 ^2 e2 

and (2.44) reduces to 

®ql ®'ll ®Q1 ^*12 ®Q2 ^12 


( 2 . 46 ) 


®q2 = ^21 ^Q1 ^21 ^22 ®Q2 


(2.47) 


with a 2^2 ®’ 2 l (2.'45)» and 


^11 = ^ ^"^dl ^ql^/^ll ^ ^22 




(2.48) 


and are now the driving point and transfer 

s s 

reactances respectively. Let 6^ and 6^ he the post- 
faxilt stable values of 6^^ and 6^ respectively* 

Define 

n = h - 

X 2 = 62 - &i 


= ij = 62 


(2,49) 


Further, let o = ~ 


( 2 , 50 ) 
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Define now the nonlinearity 


£{a) = sin(a + 5 ^ - 


(2.51) 


where = 6^ - 6® • It is to be noted that and 

4^ (2.5l) are obtained from (2.47) which shows that 
these voltages are functions of 6^^ hence cr. The 
expression for f(a) is thus seen to be much more 
complex than for the single-machine system. However, 
it can be verified that f(o) does satisfy the sector 
conditions a.round the origin, so that Popov's criterion 
can be applied in a region aroimd the origin. With 
equations (2.49) and (2.5l)» the system dynamics (2.42) 
can be cast into the form; 


Xp = x^ 


Xp = x^ 


= [-DpX^ - f(o)]/Mp 

^4 " ^ 


:2^52) 


with a as in (2.50). The system (2,52) can be transformed 
into the standard form (1.8) by the change of variables 

= Dp Xp + Mp x^ 


^2 ~ ^2 ^2 ^2 ^4 


(2.53) 
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Replacing x^-and x^ with, these new variables, the system 
can be represented in the form 

X = F X - g f(o) 

£ = -d f(cr) 


T T f 

a = h X 4- £ 2 


(2.54) 


where 


! X- 


X = ! 


i 


1 

2 


(2.55) 


and F = 


0 . 


d = 


'll 

i 

1 

f h = 

--Mi/Bi ■■ 

'-1 ^ 




g = 


1/M^ 

-I/M2 


l/Dp 

-I/D2 


(2.56) 


The transfer fimction, of the linear part of the above 
system is given by 


W(s) =.h^(sl - F)”-^ g + P^' d/s 


•"1 


.T 


(2.57) 


lefine now 

Z(s) = (2o:£^ d + P s) ¥(s) (2.58) 

Fote that unlike in Section 2.2.1, here £ is a vector. ' 

T 

£ d IS a positive scalar though, and we can conveniently 

T 

choose a such that 2a p -d = 1. As in Section 2.2.2, 
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a frequency criterion given "by (2.20) can be applied here 
also for the existence of a liapunov function of the form 
(2.17). The construction procedure then follows the same 
way as, in that section. 


2.3*2 liapunov function: 

As discussed in Chapter I, in a multimachine system 
two distinct cases, 'nOnuniform* and ‘uniform' damping, 
arissi Both these cases will be dealt with for the two- 
machihe system^*^. 

Case (i) Ilonuniform Damping: 

Let D^/M^ = and ^ ^2 ^ ^ 2 ^2. 59) 


The transfer function derived from (2.57) using the 
appropriate matrices and vectors from (2.5^), is seen on 
simplification, to be 


¥(s) = 


(1/M^ + 1/M2)s + (\M2 ^2^“l^ 

s(s+^J^) (s + ^ 2 ^ 


(2.60) 


Define 


1/M^ + iyM 2 = a ; \/E^ + X^/M^ = b 


Choose 


a 


Vs 




T 

so that 2a P d = 


1. ¥ith this, we have 


Z(s) 


b+(a+pb)s + aPs^ 

s(s"+x^5rs'"T^ 


( 2 . 61 ) 


(2 . 6l) 


( 2 . 62 ) 
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Using (.2.62), we get 

i[Z(jW) + Z(-3“)] = [X^X2(a+U,3)-b(X^+X2) + w2(ap(X^+X2^ “■ 

(a+l3p))]/[(w^+X^)(w"+Xp] (2.63) 


Therefore, to satisfy the Popov’s frequency condition 
(2.20) or equivalently for (2.63) to he nonnegative , it 
is necessary that 


^1 ^2^^^^^^ - b(X^ + X^) > 0 

(2.64) 

and 


ap (X ^ + X 2 ) - (a + bp ) > 0 

(2.65) 


Substituting for a and b from (2.61) the condition (2.64) 
can be seen to reduce to 

- (V \ ^ ° (2-P6) 

This can always be satisfied by appropriately choosing p. 
Clearly, different values of p will result in different 
V-functions. A possible choice is 

p = 1/X^ + l/X^ (2.67) 

With this choice of p, the condition (2.65) will simplify 
to 2 + ^1^2^ “ ^ which is always true. 

Indeed the left hand side of (2.65) is strictly positive. 
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Construction of V-functlon: 

Steps involved in constructing T»-function, are 
the same as in 2,2,2 except that, in the steps (1) and 
(ii), p_ d is to he used instead of mere p. 

(i) With Z(s) as given In (2.62) it can he shown that 


1[Z(b)+Z(~s)]= 


k, - s' 


(-s^+X^) (-s^+Ap 


( 2 . 68 ) 


where 


= D^D2(M^+M2)/(M^M2)^ k2 = 

(2.69) 

factorizing (2.68) in the form given in (2. 21), we have 


T(s) = (Vk-, + fk^ s)/[ (s+\ ) (s+X«) ] 


(2.70) 


(ii) Similar to (2.22), now we have 


Y = p d + h^ ■ g) 


which on substitution of appropriate vectors results in 


Y = 0 


(2.71) 


(iii) Using (2.23), li is to. he determined. Let u he 
of the form 


u = 


u. 


11 


u 


21 


(2.72) 
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Then equation (2,23)iJ.pon substitution yields 


X - ■»_, u, 

C ^ 


T(s) = 


^ ^11 ^M2 


(s + X^)(s + .X 2 ), 


(2.75) 


V/ith T(s) as given in (2,70), we obtain from (2.75) "by 
equaling coefficients of like terms 


-'^11 2/1 + '"21 ^ \ l ^2 = frl 


(2.74) 


Solving simultaneously the two equations in (2.74), we 
have 


Uii = - \ - ^2) 


U21 = " ^2^ 


(2.75) 


(iv) let the matrix 1 in the Liapunov matrix equation 
(2.24) be of the form 


^'11 ^12 


12 22 


(2.76) 


Substitution of P, u and 1 in (2.24) results in the 
identities; 

-2 = -uq 


-2 il X 
22 2 


"U 


21 


^ 12 ^^ 1’‘’^2 ^ '^ 11^21 


(2.77) 
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SulDstitution of ^ and from (2.75) in. (2.77) > will 
give on simplification: 

2X,(x^-xp = 


M;(k^ 1^, - gl, If(kiV)) 

2X2(X^ -xp^ 


M, 1, 

*12 “ (X +X )(X _X ^ 

^ ^ ^ ^ (2.78) 

1 can be verified to be positive definite. The Liapunov 
function is now given by 


V(x, a) - (£^^ + 2-2"^ ■''^^^12~ ^ 


ao^ + “ax_-^o'x«+|3 / f(a)dcr (2.79) 

1 ^ 2 ^ o , . 

where a, p and f(a) are as defined in (2.61), (2.67) and 
( 2 . 51 ) respectively. 

Case (ii) Uniform Damping: 

Consider the case when 


h/“l = ®2/“2 = 1 


( 2 . 80 ) 


In .this case it is possible to combine the two equations 

in (2.42) into a single differential equa,t j on dff.thh / 

1 T . jy OENTRal LJuti ■■■■ f\Y 

following form: 

Acc. No. /. 293 . 9 i 
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6i2 + ^ ^12 


where 


M 




( 2 . 81 ) 


= M^E^/iE^+E^) ( 2 . 82 ) 

Other symbols have the same meanings as in the previous 
case. 


Define the state variables 


^1 ^12 " ^12 

^2 = (2.83) 

Define also 

E E 

f(a) = |^[ -M^slnia + 6 ^ 2 ) - \i] (2.84) 

where 

a = x^ (2.85) 

The voltages and Eq 2 axe to be determined as before 
from (2.47). With the above variables, the system 
equations will now be: 

il = X2 

^2 = - ^X2 -- f(cy) (2.86) 

with o = x^. Equations (2. 86 )a»e exactly similar to 

the system equations (2.10) for the single machine -• 
infinite bus system except for the difference in 
definitions of the state variables and nonlinearity* 
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Hence with, a = ^ /2 and p = l/^ ^ we have a Liapunov 
function 

X 1 ^ 

V(x^,cr) = -(c + X i / f(a)dc (2.87) 

with X 2 and a as defined in (2.83) and (2.85) respectively. 
In the case when machine 2 corresponds to an infinite bus, 
(2.87) and (2.35) are identical. ¥ith the Liapunov function 
available, the stability study of the system can be carried 
out incliiding the determination of the critical clearing 
time. 


2.4 EXTENSION TO k-MACHINE SYSTEM 

The method discussed in the previous sections can 
be extended in principle to a k-machine system, including 
transient saliency. The form of the set of differential 
equations describing the system will be 

Ml S'l + D. 61 - P^l - P^l (1.1, 2,..., M) (2.88) 


is the electrical power output of the ith machine and 
in this case, the determination of this quantity involves 
considerable labour. It can be shown^"^’^ that the 
electrical power output of the ith machine is given by 

k 


ei 


Eq. 

a + I 
^11 j.i 


jVi 


E E 

sin(6 - a ) (2.89) 


where E and E^ . are the voltages back of the quadrature 
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axis synchronous reactance of the ith -and jth machines 
respectively. and are, as hefore, the driving 

point and transfer impedances respectively, and 

have the same meaning as hefore. It is to he noted that 
^ij ” “^ij - voltages Eq^, Eq 2 »»»*» 

^Qk determined from the following set of 

k simultaneous equations: 


S' 

qi 


(x' .-X . ) 


J -i- Ij 


( i— 1 ,2,...,k) 
( 2 . 90 ) 


If we neglect now the resistances as in the previous 
cases, (2.89) and (2.90) reduce to 


k -ri -ri 

P j = I — ^ sin 6. . 
ei 3=1 


(i = 1,2, ... ,k) (2i9l) 


e:. = [1 + ]V - L ^ 

n j— -1- ly 


qi 


Qi 

(i=l, 2 , . . . ,k) 


13 


(2.92). 


Thus E . (i = l,2,...,k) can he expressed in terms of 
the known voltages (i = l,2,...,k) which are' 

proportional to the field flux linkages and substituting 


for these voltages, the power P . (i = l,2,...,k) can be 

63L 
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determined for solving the differential equations (2.88). 
It is to "be noted however, that (2.92) is a set of 
k algebraic equations in (i - l,2,...,k) whose 

coefficients are functions of 6. Furthermore, when 

X J 

the order k is high, an explicit expression for Eq^ in 

terms of the known voltages S’. (i=l,2,..,,k) is 

qi 

highly complicated and when substituted into the power 
expression (2.91), results in a nonlinearity that is 
very difficult - if not impossible - to handle. There- 
fore, at the present time, no attempt is made to develop 
a liapunov function for multimachine system with saliency 
effect included. 


2.5 COECLUSIOE 

In this chapter, the saliency effect has been 
included in the stability analysis of power systems, 
in some depth. A single machine infinite bus has been 
discussed in great detail. A liapunov fimction V has 
been derived and this 7-function has been utilized 
for a parameter analysis of a typical system. This 
example illustrates clearly, the power of the Liapunov 
approach in making large number of studies with very 
little effort and time compared to the conventional 
step-by-step method. A two-machine system is then 
formulated. Both the 'nonuniform' and 'uniform' damping 
cases have been dealt with and appropriate V-functions 
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derived. Studies similar to those conducted for 
the single machine- infinite hus system can he easily 
done for this system also. The difficulty of extending, 
to a k-machine system, where k is large, is pointed 
out . 

. , The models for which V-functions have heen 

successfully developed, can be allowed to' include 
further details such as flux decay effect, voltage 
regulator and governor dynamics, ¥e will see the 
case of a single machine- infinite hus system with 
such details in the next cha.pt_er‘ 



CHAPTER III 


SINGIE-MACHINB SYSTEM WITH IJyjPROVlD MACHINE MODBI 

3.1 lETRODUCTIOH 

In this chapter a system with a single machine 

connected to an infinite bus, referred to hereafter as 

a 'single machine system' for brevity, is investigated 

further. As mentioned at the end of the last chapter, 

the synchronous machine is modelled to include further 

details such as flux decay, voltage regulator and 

governor dynamics. At first a model incorporating only 

flux decay and voltage regulator is studied. It is 

observed that the system now has multiple nonlinearities. 

furthermore, the nonlinearities belong to the class of 

'multiargument' type because of which a direct application 

of the matrix version of Kalman's construction proce- 
38,39 

dure is not possible. However, a modification of 

this construction procedure by incorporating the results 

40 

of Desoer and Viu provides us with a method to obtain a 
Y-function. Using this V-function, a parameter analysis 
is carried out to investigate the effect of the main 
parameters of the voltage regulator. The same 
construction procedure is then extended to a system 
model wherein saliency effect and governor dynamics are 
also included in addition to the flux decay and voltage 
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regulator action, liie Liapunov function derived for this 
model is seen to he most general and naturally the Liapunov 
functions derived earlier happen to be special cases of 
this V-function. 

5.2 SINGLE ILiCHINB SYSTEM ^wITH ELUX DECAY AND VOLTAGE 
REGULATOR 

Several authors^' have dealt with 
the effects of flux decay and voltage regulator action on 
the transient stability of synchronous generator. Indeed 
these effects are being included in many of the digital 
simulation studies. It is only appropriate therefore that 
these effects are included in stability studies via 
Liapunov method. 

The dynamic equation governing the flux decay is 
simple and straightforward. Inclusion of this in the 
system equations therefore poses no particular problem. 
However, the nonlinearities involved in the system 
dynamics now get changed to the ’multiargument' type, 
the implication of which has already been mentioned in 
Section 5*1« 

The voltage regulator action on the other hand, can 
be represented in several ways. The exact model of a 
voltage regulator renders the problem difficult for 
solution through Liapunov approach. We will therefore 
approximate it by a model having an exponential type of 
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response . The finite time that elapses before the 

voltage regulator begins to act following a fault, is 

neglected partly because it is very small and partly 

for the sake of mathematical simplicity. With this 

representation of voltage regulator, the system equations 

become nonautonomous . Siddiqee used this model as it is 

and his attempts to derive a V-function through known 

12 

methods of construction seems to have failed . Indeed, 
it appears that the well-established construction pro- 
cedures that have been successfully employed to derive 
T-functions for autonomous systems, cannot be extended to 
nonautonomous systems. Fortunately however, with the 
above .model for voltage regulator, it is possible to 
convert the nonautonomous problem of the power system 
into an autonomous one, by a change of variable. This 
enables to derive a Liapunov function using a modified 
version of Kalman's construction procedure. 


3.2.1 Formulation of the Problem: 


Neglecting saliency and resistances the rotor 
dynamics of a single machine connected to an infinite 
bus can be written as 


M 6 + D 6 


J- 


B 


'm 


x^+xl 
e d 


sin 6 


(5.1) 


where is now the infinite bus voltage. The flux 
decay effect can be represented through the variation of 
which is given by 
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E’ = [E ^ (E’ + (x,-xj )i, 
q “-ex q d d^d-'^o 


(3.2) 


where 


Suhstituting (3'*3) into (3*2) and simplifying. 


B 


E' 

1 


ez 


% ■ ^e+^d (=^e«P % 


COS 6 


(3.4) 


The voltage regulator action on the exciter voltage is 

54 

approximated by 


E 

( 

where 

E 


= E - Z E e 
ex c VO 


Eo(l + K^) 


“E4t 


(5.5) 


(3.6) 


with E^ as the initial exciter voltage at the time of 
fault. Define 

= 6 - 6'^ 

« • 

X 2 = = 6 


z, = E’ 


^4 " ^ex ' ® 




(3.7) 


where e is the post fault steady state value of E^. 


Define also 
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X +x, 
e d 


\ + X' Ti' ' 




T]^ = 1/T^ ; Ti^ = 1/T^ ; = Eg/(Xg+x^) (5.8) 

Using (3*7) and {3-8)» eqiiation (3 *4) can Ue rewritten 
after some manipiilations as 


X 


^ + T]^x^-r) 2 (cos 6^ ~ cos(x^ + 6®)) 


from (3*7) 


^4 = “"U ""4 


( 5 . 9 ) 


(5.10) 


Define now 

K ' P 

(<^2 sin(a^ + 6®) - 

= TI 2 CCOS 6® ~ cos(a^ + 6®)] 

where 




( 3 . 11 ) 


(3.12) 


Making use of the above relationships, the overall system 
can he represented hy a set of first order differential 
equations as ; 

±1 = x^ 


= -i X^ - 

% = "^ 1^3 + ^ 3 \ ^ 


^4 ~ ”''^4^4 


( 3 . 13 ) 
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Defining a new state variable ? as 


5 = X 4- 


(5-14) 


and by a change of variable, the system can be represen- 
ted by the alternate form 


- . i 


X 


2 ! 


1 

« 1 


5 


_ 

i 


-X 0 

0 -T) 


1 


0 


0 0 ^T] 



^2 


1 - 

o 

H 

L 


^3 

- 

0 1 


''4 


0 0 

j 



) _ 


?=-[! 0] 


fl(2) 

f^Cai) 


-M/D 0 0 

0 10 


fl(£) 

f2(^l) 



^2 


M/D~ 


x^ 




^4i 


0 


(3.15) 


This is now in the standard form 


y = F y - & f( 0 ) 
k = -d^ f (o) 

■ 0 = y + p C (3-16) 

The nonlinearity f (£) above has one of its components 
f^(p) having two arguments. Furthermore, fj_(o) violates 
the sector condition away from the origin. Nevertheless, 
it is possible to develop a Liapunov functian which 
satisfies the required sign definite properties in a 



63 


sufficiently large region around the origin as has been 
established by the numerical example given in Section 3*4. 
The transfer function matrix for (3«16) is given by 

¥(s) = H^(sl - i “ P (5-17) 

s 


Comparing (3.16) and (3.15 )» and substituting the 
appropriate matrices and vectors in (3.17), the transfer 
function matrix for the system (3 .15) can be derived to 
be 

0 

S (s+ ^'- ) 


V7(s) = 


(3.18) 


0 


1 

(s + h^) 


It is to be noted that system (3.15) is a multinonlinear 
system with two nonlinearities and has a transfer 
function that is a matrix instead of a scalar. Thus 
a matrix version of Kalman's construction procedtire has 
to be used. But a direct application of even this 
procedure is difficult for this case, as the nonlineari- 
ties here belong to the class of ' mult iar gum ent’ type 
nonlinearities. The difficulty is overcome by suitably 
adapting Kalman's construction procedure by incorporating 
some aspects of the work of Besoer and ^ 
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3.2.2 Liapunov Function for Systems with. Multiargument 

55 

ITonlinearities : 

Consider a nth order dynamical system belonging to 
the simplest particular case characterized by (3.16), 
where 'jy is a (n-l) state vector, C is the nth state 
variable, F is a strictly Hurwitz (n-l )x (n-l) real 
matrix, G and H are (n-l) x m real matrices, d and jp 
are m-vectors. The m-vector fCn) belongs to a class 
of multiargument nonlinearities, which implies that 
each of its components f^, may be a fiinction 

of all or some of the m arguments 

Now, the following assumptions are made as in Lesoer 
and 

al: f(£) is continuous and maps into R^ 

a2s For some constant real symmetric matrix R and 

T 

such that R _pd is symmetric 

rp rp -m 

2 R^d^ f(o) > 0 for all a e R 
and 

f(a) = 0 ■ if o = 0 (3.19) 

1 ^ 

a3s There is a function e C mapping R^ into R 
such that 10 for all o e R°^ , with 

Y^(0) = 0 

and for some constant real matrix Q 

for all real o e r“^ 


( 3 . 20 ) 
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Now define 

Z(s) = (2RP_ + Qs)¥(s) ( 5 . 21 ) 

A theorem, which is a modified version of the one due to 
47 

Anandamohan is given "below and proved in the Appendix B. 
Theorem ; 

Bor the system ( 5 .I 6 ), let the assumptions al-a3 
"be valid, hith constant matrices E and Q defined as in 
a2 and a3 and Z(s) as in (3.21), the frequency condition 

4-[Z(jw ) + w)] > 0 for all realw (3.22) 

is necessary and sufficient for the existence of a 
Liapunov function 

V = L y + (o - y)^ R(o - y ) + v^(a) (5-23) 

with the derivative a.long the solutions of the system 
equations (3*16) "being given "by 

V = -I y + r f(o) jl ^ - 20^ Rp d^ f(R) . (3-24) 

where 1 is a (n-l)x(n-l) symmetric positive definite 
matrix, "U is a (n-l)x m matrix and T is a m x m matrix. 

Procedure for Construction of "^-function: 

The construction of Y-function for the system 

( 3 . 16 ) essentially involves choosing the matrices R and Q 

first, and then arriving at the matrix 1 through the 

34 

systematic procedure due to Kalman . Although there is 
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no precise rule for clioosing R and a close examination 
of the formulation given above gives some general guide- 
lines to arrive at these matrices. For instance, Q and 
are related through equation (3.20), so that, once a 
suitable Y^ is chosen, so as to satisfy assumption a3, 

Q is fixed. As in the classical Popov- type Liapunov 
function, Y^ can be chosen to be of the general form of 
'an integral of nonlinearities’ type. But because of 
the conditions imposed by a3» it may not always be 
possible to include all the nOnlinearities of the problem, 
in Y^. It may be possible to have several Y^»s and 
therefore, choice of Y^ is not unique. In such a case one 
would get correspondingly many Q’s and hence many Y- 
functions. It would appear that inclusion of more number 
of nonlinearities in Y^ would yield better Y-functions, 
implying thereby larger stability regions, as in such 
Y-functions the system would have been reflected to a 
greater degree of details. Similarly, R can, in general, 
be chosen arbitrarily, but must be such that it should 
satisfy assumption a2 and together with the chosen Q, it 
also should satisfy the frequency condition (3.22). Once 
R and Q matrices are chosen, the rest of the construction 
procedure is based on the proof of the matrix version of 

38 

Ealman-Tacubovich Lemma discussed by Farendra and G-oldwyn 

39 

and Narendra and Neuman and used earlier by Pai and 
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Anandamohan^^ and Pai^^ for multimacliine systems. 
Essentially, tlie following steps are involved: 

(i) The matrix M(s) is solved for, from the relationship 

i[Z(s)+Z^(-s)] = M^(-s)M{s) (3.-25) 

This ca,n he clone, in general, by spectral 
factorization^^”^^. 

(ii) The matrix F is determined from 

iQ(H^G +£d^) + + d jP^)Q^ = r^r (3.26) 

(iii) The identity 

M(s) - r = -U^(sl - P)"^ G (3.27) 

is solved for the elements of the (n-l) x m 
matrix U'. 

(iv) Einally the Liapunov matrix equation 

P^L + L P = -U TJ^ (3.28) 

is solved for the (n-l) x (n-l) matrix 1. With 1, 

R and V^(£) known, the Liapunov function can be 
written as in (3.25), with its derivative given by 

(3.24). 

3 . 2.5 Application to the Problem of Section 3.2.1: 

In the problem discussed in 3.2.1, there are^ 
two nonlinearities f 2_(o) and It^eally one would 

wish to include both these nonlinearities in Y-^{o) , But 
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it is Observed that this leads to difficulties. The 
only convenient choice in this case turns out to be 

\ = q. / f^Co) d cr 3 _ (3.29) 

with q. > 0. Now applying the ’gradient condition’ 
( 3 . 20 ), the matriz Q is derived as 


Q 


r 

j q. 

I 

1 ° 

L 


qK. 


1 




R can be in general of the form. 


(3.50) 


"'ll 

^12 

^12 

^22 


m 

Comparing (3.15) with ( 3 .I 6 ), the matrices F, G, H 

and the vectors P. , d can be identified for the system 

under study. It is observed that with the £ and d 

T 

of the system, to make R£d a symmetric matriz, R 
should be necessarily a diagonal matriz. Thus with 
the off-diagonal elements of R assumed zeros, the 
condition (3.19) requires that r^^ f;i_(2) -2. 0, i.e. 

riiai[Ki(a2+e) sin(a^ + 6®) - Pj/M > 0 (3-32) 

For the problem under consideration, the above condition 
may get violated in some region around the origin (for 
ezample, for a fized positive value of a^) , which 
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in t-urn would adversely affect the sign definite property 

♦ 

of V. To avoid such a situation, R is chosen as 261*0 
matrix j 

R = 0 (3.53) 

R and Q have thus "been decided upon. With Z(s) as 
defined in (3.21), it can he shown that 


ilZii “ ) + Z^(-D w )] = 


(3.34) 

'^2 "being positive. Row using (3.25), we 
get on factorization 


1 X 


2 . 2 

CO + X 


0 


0 


((lZ^/Mr]2)«^ 

2 2 
“ + ri£ 


> 0 


M(s) 


f(q X )/(s+ X ) 


0 


sf ( qE^/Mh^ ) / (s+R^ ) 


( 3 . 35 ) 


Next solving (3.26), the matrix V is obtained as 

1 


r = 


0 


0 


0 l/'(qK^/MTl2) 


(3.36) 


let the 3^:2 matrix U he of the form 
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^11 

^12 

^21 

^22 

U31 

^32 


( 5 . 37 ) 


Substituting the appropriate matrices derived above into 
(3.27) and solving for U, we have 


^11 ~ ' "^12 ^21 ° ' 


U 22 = t,3_ /{,Kj_/Mr,2) 


(3.38) 


The remaining elements u^^ and u^^ arbitrarily 

chosen. Choose 


u^l = 0 ; u^2 = - f(q%/MP2) 


(3.59) 


(Note; This choice of u,t and u^^ has been made 
deliberately, so as to obtain a diagonal 1 matrix, 
for simplicity. In general however, 1 need not be 
diagonal). 


¥ith the elements of U given as in (3*38) and 
(3*59), the Liapunov matrix equation (3*28) is now solved 
for and the matrix 1 is found to be 


1 


q/2 0 0 

0 qK^p^/CaMn^) 0 

0 0 qTL^r]^/ {8m ^r)^) 


(3.40) 


which is clearly positive definite. 
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Nov V^(£) is calculated from (5*29) after sulDStitution . 
for- f^(£) from (3.11) together with e sin 6®, 

and is found to be 

K. o o 

= <1 [|;^( ff^+e) (cos 5 -cos (cr^+6® ) )- ^ ® 

(3.41) 

Without loss of generality, q can be chosen as unity. The 
Liap-unov function will now be given by 


7(x, £) = i ^2 + 


^1^1 
2Mp^ 3 


^3 ^1 


^4 + 73_(o;) (3.42) 

r 

X,-. 


8Mp2'n4 

where given by (3. 41) with q = 1 and x = 


The time derivative of V along the solution of the system 
equations is seen to be 


K. 


V = 


T)^- K, 


1 ( —■ \2 h i. 3 


1 E-iT]- 

1 2 1 3 * 

- X. + x^x 




3 4 


(3.43) 


which is seen to be negative semidefinite in the region of 
interest. 


The above Liapunov function V is seen to be 

12 

different from that given by Siddiqee and is foimd to 
give more satisfactory resiiLts for the numerical example 
given in Section 5.4. In fact, Siddiqee 's V-function was 
observed to violate the required sign definite properties 



for certain initial operating points, even in a region 
very close to the origin, thus making it unsuitable for 
certain operating conditions. On the other hand, the 
V-function derived above, remains valid for all likely 
initial operating conditions, in a sufficiently large 
region, which contains the region of interest. 

3 . 2.4 Equilibrium Points and Stability Region: 

It an equilibrium point of the system under 

• ♦ • pg 

discussion, = X 2 = = 0. It can be shown that 

the solution of the following equations determine the 
equilibrium points of the system: 

. K.ri 

^m " ®ex ^ ^ 26 = 0 (3.44) 

= (^3 ®ex ^2 ^^/^l (3.45) 

Solving these nonlinear algebraic equations, the stable 
prefault equilibrium point is calculated with the 
corresponding system parameters. Similarly the post 
fault stable and unstable equilibrium points are computed 
with the parameters corresponding to the post fault 
condition. The stability region is found from the 
inequality 

V < c (3.46) 

where C is the value of V evaluated at the unstable post 
fault equilibrium point. The problem above being that of 
a single machine, the identification of the stable and 
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unstable eq_nilibriuin points is straightforward. The 
critical clearing time is computed by numerically 
integrating the system differential equations (5*15) 
for the faulted condition, starting from the prefault 
operating point until the solution trajectory reaches 
the stability boundary G. 


3 . 2.5 Augmented V-function: 

The Liapunov function (3 ‘42), can be made more 

general by augmenting the quadratic terms as done by 
27 

Mans our . Thus we can have 


Y = 


+ ix; 


2Mri, 


X^ + 

3 




x^ .+ 




(5.47) 


The choice of a and p have to be such that the sign 
definite properties of V and V are not violated in the 
region of interest. It can be verified that the 
quadratic form in Y will remain nonnegative if a ^ 6^/2. 
Additional restrictions have to, be imposed by examining. Y. 
It is observed that for this problem, it is not possible 
to make Y negative semidefinite in the whole state space 
with these additional terms present. The objective of 
introducing these terms is, however, clear. It is to 
obtain a larger region of stability which is closer to 
the actual stability region. Therefore, the best one can 
aim is to select a pair of a and p which will achieve 
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tills objective to the extent, possible, without violating 
the sign definite properties of V and Y. This calls for 
certain amount of search as evidenced by the numerical 
example given in Section 3*4. 


30 G-ENERAIIZED MODEL OP |IMG1IHMACHINE SYSTEM 

In this section we will consider a single machine- 

infinite bus system with the synchronous machine modelled 

to include saliency, flux decay, voltage regulator and 

governor dynamics. Clearly, we should arrive at a Liapunov 

function that is most general, of which the previous 

V-functions must be mere special cases. Indeed it will be 

shown to be so. Although excellent models with great 

details are available in the literature for governor 

representation^^* a speed governor with one time 

20 27 

constant will be considered here for simplicity * 

Higher order representation does not really pose any 
problem. 


3.3*1 Eormulation of the Problem: 

The governor action can be approximated by an 
eq^uivalent first-order system as follows; 


T„ + P^ = <5 

g m m mo g 


(3*48) 


where is the instantaneous mechanical power input, and 
P is the post fault steady state value of P . K is the 
governor amplification factor and T is the equivalent 
time constant of the governor sysleni'. 
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Define now a. state variable 

x_ = P - p 
5 m mo 


(3-49) 


With the state variables to as in Section 3* 2.1, 
the dynamic equations of a single-machine system with 
saliency, flux decay, voltage regulator and governor 
dynamics can be written as 


^1 = 

^2 



X2 = 

- X x^ 

+ M ^5 “ 

fl(s) 

^3 = 

-^3 

“f" "" 






^5 = 

“g 

5 


g = 

K /T 
g^ g 

and a 

= 1/Tg. 


(3.50) 


the nonlinearity f^ (o) is now given by 


f^(c) = [K^(a2+e)sin(o^+6® )-K2Sin2(cr-, +6®)-P^]/M 


m' 


(3.51) 


where 


X, 




; P = K-, e sin 6® - K, sin 26® 
^ m 1 ■ 2 


(3.52) 

All other symbols have the same meanings as in 
Section 3.2.1. System (3-50) can now be transformed 
into the standard form (3.16) by a change of variable 
aD+g , x^ 


20 


(■ 


Ma ^^1 


)x-| + X, + 


Ma 


(3.53) 



76 


which results in 





(3.54) 


For the above system, the transfer function matrix for 
.the linear part obtained by substituting the appropriate 
matrices and vectors from (3-54) into (3.17) > is seen 


to be 


W(s) = 


ah+g aD+g 

s ( s+ ^) ( s + a ) 

0 


0 

s^+(a+x)s + (aD+g)/M 
(s+>^) (s+p^) (s+a) 

(3.55) 
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Note tha-t the above transfer fimction reduces to the one 
derived in Section 3 •2.1 when governor action and saliency 
are neglected. 


3.3*2 Derivation of Liapunov Function V: 

Following the procedure outlined in 3*2.2 a 
V-function for the above problem will be derived; As in 
Section 3.2.3» let 

\ = q. / (3.56.) 

0 


.A. 

where f^Cs) is ixow given by (3.51). On integration we 
get as 


"7^(2) = ^^(cr^+e ) (cos6®-cos (a^+6® ) )- •~(cos26®-eos2 (a^+6® ) ) 


P 

^ ^ 1 
M 1-* 


(3.57) 


With V^(0) as given by (3.57); the application of the 
gradient condition (3.20) yields the matrix Q as 


Q 


0 


0 


qK3_/(Mp2) 


(3.58) 


Interestingly enough, the matrix Q turns out to be the 
same as the one in Section 3.2.3. Again for the same 
reasons discussed in 3.2.3» R is chosen as a zero matrix 
for this case also. Now with Z(s) defined as in 
(3.21) and W(s) and Q given by (3*55) and (3*58) 
respectively, we haves 
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MZ(jco)+Z^(-^u3)] = 


^110 

(a) ^+A ^ ) (u ) 


0 


— [ 

Mt] ^ 


0 


(0 ^ J.]r w'^+k CO' 

^^224 ^^222 


2 ( a32+X2) (a)2_^^2) 


(3.59) 


where 

k 


XaD-g(aHA ) 

= ^ ; k 

aD+g 


a^ D X 


aB+g ’ 


^224 ^ ^“g/M ; ^222 ^ X^+|( XT}^+ar]^+ax ) 

(3.60) 

So as to sa,tisfy the frequency condition (3.22), we shall 
impose the following two conditions: 

X a D - g(a + X)_> 0 

a2 +X 2 _ g/M 1 0 (3.61) 

How the rest of the construction procedure follows 
on the lines indicated through (3*25) - (3*28). Factorising 
•|-[Z(s) + Z^(~s)] as given hy (3.25)» we get M(s) which is 
seen to he 

— ■ 


M(s) = 


0 


m22(s) 
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where 


m22^®) = f(K]^q/MTi2) s(s+ V'[-|^ - V(' -— — ^)])(s + 

^[_|24 ^ ^-(.,. 2 , 24. ^ . ^ . ^ )])/[ (s+X) (s+n^) (s+a)] 

( 3 ^ 62 ) 


Now with M(s) as given "by ( 3 . 62 ), making use of the 
relation (3-26), the matrix f is obtained as 

r 0 0 1 


r = 


(3.63) 


0 


\fi- 




Mr). 


Assume the 4x2 matrix U to be of the form 


Uii 

U12 

U21 

^22 

^31 

U32 

H 

'cT 

P 5 

1 

^42 


(3.64) 


With this form of U and substituting the appropriate 
matrices into the identity (3.27)? the elements of U are 
solved for and we get explicit values for the following 


elements as 
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: 


"^12 = ^21 " "^42 ° ^ 


"^41 Vq. ( T^2.12^^® ^ 


Z^q 


'^22 = ''(mti 


■) 


aX T] 


1 


2 (a X+ g/M) 


(3.65) 


The remaining elements "^32 arbitrarily 

chosen, depending upon which choice, different I matrix 
and hence different V-fimctions will be obtained. Por 
the sake of simplicity however, we will so choose these 
two elements, as to obtain a diagonal I matrix. This 
requires the choice of and u^^ 


u 


31 


0 


u 


32 


P, K. q 
_2 ][( — i_) ( 
2 


aX 


+g/M 


) 


( 3 .. 66 ) 


Using the above U as given by (3.65) and (3.66) and 
solving the liapunov matrix equation (3-28), the 
L matrix is obtained as 

L = Diag(^u^) (i=l,2,3,4) 

w^here 


11 


22 


'33 


= q[x aI)-g(a+X)]/[2 X(aD+g)] 


K q p. 


2Mp, 




a X 2 
] 


■a'x+g/M 

a X 


aX+g/M 


) 


q{fk 




110 


afk^io) 


-^44 


2g^a 


(3.67) 
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where ^22.0 ^112 given in (3.60). The liapimov 

fiinction for the power system under consideration is given 

■by 


v(£,2) = hl^2 + * '■44^5 * h<2) (3-68) 


With 1 given by (3.6?) and as in (3.57). It is 

easily verified that this Liapunov function reduces to 
the Y-function derived in 3*2.4 when the governor dynamics 
and saliency are neglected. On the other hand, if we 
consider only saliency and neglect flux decay, voltage, 
regiulator and governor dynamics, the V-f unction given 
above reduces to the same as the one derived by 
Mansour (except for a constant multiplier) by the 
energy metric algorithm. The Y-function for this case 
as derived from (3.68) differs slightly from the one 
derived in Section 2.2.3. This can be attributed to the 
difference in the choice of R and L matrices in this 
section and the corresponding quantities p and ^ in 
Section 2.2.3. Nevertheless, this is a perfectly valid 
Y-f unction. In general therefore, the Y-function given 
in (3.68) is a very general one and if we wish to neglect 
any of the details that have been included in the above 
model, we merely set the corresponding state variable 
as zero and make appropriate changes, in the relevant 
coefficients. As in Section 3.2,5, the Y-function can 
be made more general by augmenting the quadratic form. 
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3.4 1 JnJMERICAI EXlD-rPIE - MODEL «ITH ELUX DECAY AID 
VOLTAGE REGULATOR 

Consider the power system shown in Figure 2.4, 
hut with the following specifications: 

H = 2.76 M.J/TIVA ; = 1.15 ; x^ = 0.3 ; 

The generator*, rating : 25 MVA ; 60 Hz. 

Each of the double circuit transmission line has a 
reactance of 0.2 p.u. The prefault excitation voltage 
is taken as 1.22. T' = 6.6 sec. A symmetrical 

oA 0 

3-phase fault is assTjimed at the middle of one of the 
parallel transm.ission lines. 

Part ( i ) : 

With the aid of the Liapunov function (3.4-2), 
the above system is investigated for its transient 
stability properties with the synchronous machine 
modelled in the following three ways; 

Case (a): Flux decay and voltage regulator action 

are considered with T^ = 5 Sec. and = 1.- 
Case (b): Only flux decay is considered and the voltage 
regulator action is neglected. 

Case (c): Both flux decay and voltage regulator action 
are neglected. 

The critical clearing time for the various cases 
are estima.ted for different prefa-ult conditions (This 
is done by varying and for various damping coefficients. 

The results are given in Table 3.1* 



83 


Table 3»1» Critical Clearing Time t 
(T = 5 sec. ; K = 1.0) 

Y V 


i 

m 

D 

Criticsl 

clearing time 

t ■■ 

c 



Case(a ) 

Case {b) 

Case (c) 


0.0 

0.36 

0.2 

0.62 

CO 

• 

o 

0.002 

0.377 

0.22 

0. 66 


0.004 

0.378 

0.222 

0.72 

. - 

0.006 

0.38 

0.225 

0.8 . 

( 

I 


0.0 

I 

0.22 

0.12 

0.312 

H 

• 

o 

0.002 

1 0.234 

1 

0.131 

0,.32 


0.004 

0.237 

0.132 

0.325 


0.006 

0.24 

0.133 

0.331 


1 

0.0 

0.14 

0.04 

0.201 

1.2 

0.002 

0.152 

0.05 

0.204 


0.004 

0.154 

0.051 

0.207 

„J 
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The effects of the voltage regulator parameters 
and on the critica.! clearing time are studied for 
Case (a). Typical results are given in Table 5.2 and 
Table 5*3 respectively. 

Part (ii): 

Kext the augmented Liapunov function given by 
( 3 . 47 ) is used to study the case (a) above. In ' 

Ta.bles 5*4 - 3*6, the quantity T indicates the time 

« 

after which the required sign definite property of V 
is violated. Table 3*4 shows the results with p assumed 
to be zero and a varying. Now keeping a at a fixed value 
some positive values of p are tried and the results 
given in Table 3.5. Finally some negative values for p 
and positive values for a are used to get the critical 
clearing time and the computed values tabulated in 
Table 3*6. 

Discussion on the Results; 

The numerical example given in Part (i) clea,rly 
illustrates the effects of flux decay and voltage 
regulator action. While flux decay is observed to reduce 
the critical clearing time, the voltage regulator action 
tends to neutralize this effect. The degree of neutraliza 
tion depends on the parameter and T^ although the 
clearing time is seen to be more sensitive to than to 
T^. Increasing the speed of response of the voltage 







Table 3.3: Effect of T on t 

V c 


(K^=l; 1=0.002) 


T 

V 

t 

c 

C.5 

0.256 

1.0 

0.2535 . 

1.5 

0.2515 

2.0 

0.249, 

2.5 

0.247 

3.0 

0.244 

3.5 

0.242 

4.0 

0.239 
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regulator does not increase the critical clearing time to 
a significant degree. 

In Part (ii) , the general augmented Liapunov 
function (5.47) is used to analyze the same system. It 
is seen that hy completing the q_uadratic form to include 
all the original state variables, it is possible to 
improve the results. The coefficients associated with 
these additional tnrms however, have to be chosen after 
some search. The coefficients cannot be fixed by having 
an eye merely on the V~function. Their effect on V 
also has to be kept in view. A large number of combi- 
nations of a and p were tried. It was observed that 
very small values of these coefficients did not affect 
the sign definite properties of Y and V significantly, 
but then the improvement achieved in the critical clearing 

time t was too small. On the other hand, very large 
c 

values of a and p were seen to render both V and Y 
indefinite in a region too close to the post fault equili- 
brium point. In the numerical example above a positive 
P did not improve the results. Instead it made the results 
werse compared to that with zero value of p . On the other 
hand, a negative value for p gives more accurate value 
of t^, provided it is chosen judiciously. Perhaps the 
best way to arrive at a proper pair of a and p would 
seem to be to first choose a value of a keeping p as 
zero, and then keeping' this value of a fixed, a search 
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is made for ( 3 . Clearly, the augmented Liapunov fimction 
will req.uire more number of computation due to the 
search process required to fix the coefficients a and p. 
Viihether the additional effort and computational time 
involved are worth or not has to he assessed in terms 
of the improved values of t obtained. Using this 
general Liapunov function also, parameter analysis can 
be carried out. 

3.5 COFCLUSION 

A systematic procedure for constructing a 
Liapunov function for a single machine connected to an 
infinite bus, with flux decay effect and voltage regulator 
action included, has been developed in the first part of 
the chapter. A more complete model of the sjmchronous 
machine to include transient saliency and governor 
dynamics also, is taken up next and a generalized 
Liapunov function derived. The ?-function derived earlier, 
is seen to be a particular case of this generalized 
Liapunov function. An augmented Liapunov function ofjtained 
by completing the quadratic form is also discussed. The 
numerical example illustrates the ease with which the 
various aspects of the transient stability problem of 
systems with greater details of modelling, can be investi- 
gated. Needless to mention that generalized Liapunov 
function can be utilized for a complete parameter analysis 
to include the parameters of the governor also. 
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MJLTIMACHIl'E SYSTEM VJITH ELEX DECAY EPEECT 

4.1 INTEODUCTIOE 

This chapter is concerned with the formulation of 
the multimachine problems to include flux decay effect, 
in a form suitable for stability analysis through the 
Liapunov approach and development of suitable 7-- functions. . 
Although multimachine systems with simple machine model 
(a constant voltage behind a transient reactance) ha,ve been 
popular in earlier studies, extension to incorporate 
machine models with greater details are desirable. It 
would indeed be desirable to include all the refinements 
diseussed in Section 3-3 f in a multimachine system. 
Unfortunately, the problem then becomes so complex, that 
it is practicallj?' impossible to handle with the techniques 
known at the present time. We will deal with a two-machine 
system first, and then a three-machine system, with flux 
decay effects considered. The results of these two cases 
would then be extended to a general k-machihe system. It 
is first of all noted that as in the ease of the single 
machine system considered in Chapter III, the nonlineari- 
ties are of the multiargument type. It is evident that 
Kalman's construction procedure cannot be used to derive 
a 7-f;mction. An attempt was made to use the modified 
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version of Kalman’s construction procedure developed in 

Chapter Illy hut without success* One therefore has to 

fall back on other techniques to generate a V-function. 

A number of techniques such as variable gradient method^^®^^ 

18 *18 

Cartwright method , Aizerman’s method , Zubov^’s 
11 26 

method * , method of integration by parts and method of 

13 

Infante and Clark have been used earlier on simpler 
situations. In the present work, a version of ’ integration, 
by parts* method^^’^^ is used to develop appropriate 
Liapunov functions. 


4.2 TWO-MACHINE SYSTEM 


formulation of the Problem: 

Consider two synchronous machines connected 
together through a transmission line. Neglecting 
resistances and transient saliency, and assuming constent 
input power and constant exciter voltage the equations of 
motion of the two-machines can be written as 


*1 



^1 ~ ^ml 

®1®2 

^12 

sin (6^-62) 

“s 

^2 


^2 "" ^m2 

- ^1=2 

H2 

sin(62''6^) 


(4.1) 


where E^ and E^ are now the voltages back of transient 
reactances for the machines 1 and 2 respectively and 
are assumed to be proportional to the field flux 
linkages. short circuit transfer susceptance 

between the internal buses representing the voltages E^ 
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and E^. The rates of changes of and E^ due to flux 
decay effect ca,n now he written ass 



(E 


ezi 




(1=1,2) 


(4.2) 


is the exciter voltage referred to the armature 

circuit and Ej^ is a voltage corresponding to field 

current for the ith machine. Working on similar lines as 

44.2 

in Crary it can he shown that 


E 


E . 


li 


[1+ (Xdi-^ ? ) /X. JE, ~ (x^ , cos6 


11 


di -di^Z.^— -ij 


(i,j = 1,2 ; d ¥ 


•■(4.3) 


where and X. . are the driving point and transfer 

J- J 

reactances with transient reactances.. On substitution 
of (4.3) into (4.2), we have: 


B. = E , /T* .-a. .E. + a. .E . cos6 . . (i, 3 

1 exi' 01 11 1 13 3 13 


where 

a 

a 


ii 

ij 




1,2; jf^l) 

( 4 . 4 ) 


14 


- 5 ) 


Let e^ and e2 he the post fault steady state values of 
and B2 respectively and 6 ® and 6 ® he the post fault 
steady state values of 6 ^ and 6 ^ respectively. It is 
easily seen that 

^ml “ ■^la *12 “ '^ma 
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where Y^2 corresponds to the post fault condition. 
Define the state variables s 



X4 = X2 

^5 = ®1"®1 

Xg = E2 - 62 (4^6) 

With (4.6) substituted, ( 4 - 4 ) on simplification, can be 
written ass 

E4 = “®'11^5 ■" ®'12^®2 "■ ^^6 ®2^ cos6^2^ 

^2 "^22^6 ~ ^21^ ®1 °°®^21 ” ^^5 ®1^ cos624^] (4.7) 

Now define 

c?! = Xf - X2 

02 = - 

0^ = Xg (4.8) 

and also 

f^(o) = (02 + ep)(o'^ + ® 2^'^12 ®4n(04 + ^12^ " “ml 

^2(0) = ^12 ” ^^3 cos(04 + 

f^(0) = 9-21^®! *^12 ” ^^2 ®1^ cos (04 + ^p2^^ 

(4.9) 
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With ( 4 * 9 )s { 4 » 8 ) and (4.6), the complete system differen- 
tial equations ( 4 «l) and ( 4 » 7 ) can he put in the form: 


^1 

^3 



m- 

X2 

= ^4 



X3 

= 

^3 

- f^(o)/M. 


= - X2 

^4 

+ fi(£)/M, 


= -a^2. 

^5 


^6 

= -a 22 

^6 

- ^3(2) 


Let us introduce a change of Tariahles now, defined by: 



- D 2 _ 

+ 

“1 

^3 

h 

= ©2 ^2 

+ 


^4 


(4.11) 


With these variables replacing and the system 
above can be transformed into the standard form; 

X = L X - G f(o) 
i = -I f (a) 

cr=H^x + PC (4# 12) 


where 


X = 

^3 

^4 

J £(£)= 

fq(£) 

11 

^1 

; 2 . = 

"2 1 


^5 




^ 2 

_ 
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F = 



0 

0 

0 



1 

H 

H 

0 

0 

0 

-h 

0 

0 



-l/M^ 

0 

0 

0 

0 

”^11 

0 

o 

f 

G = 

0 

1 

0 

0 

0 

0 

”^22 



0 

0 

1 _ 




1 0 

0 







F = 



a 







-1 0 

0 

■ f 






0 

0 


1/D^ 

-1/D2 

0 

0 

1 

0 

f p = 

0 

0 

0 

0 

0 

1 


0 

0 


(4.15) 


Liapunov Function for the Two-Machine System: 

Assume tentatively a Liapunov function of the 

form 

2 .) = + 2-22 X* + === 4 - 9 . x’ 

Oi 

+ q / f^(a) d cr^ (4.14) 

n 

where q is yet to he determined. The inclusion of 
integrals of other nonlinearities was found to give 
rise to some difficulties in ensuring the sign definite 
property of the time derivative of the Liapunov function. 
The positive definite property of can he established 
along the lines similar to the single-machine -system. 
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Novj- taking the time derivative along the system 
equations (4 .12) and on some simplification, we have 

^0 = ^ 2^22 ^ 1 ^ 4/“2 

+ 2*33X5*5+2144 xgig + a - a 

B g 

+ q_ Y^2 (cf3+e2)[cos + 6^2 )] 2;^ 

S * 

+ q. Y2_2(°'2 ■*■ ®p)[ cos6^2 ^12^^ ^6 ( 4 . 15 ) 

Choose 


\l = ^1*^/2 ; ^22 " 

Viith (4.16) substituted in (4.15) aud on some rearrange- 
ment, we have 


V 


-2 ?-T. A, - j2 S' X ^ X 


■'ll 1 "^5 
2 i 


' 2 

.' 22 2 ^4 




^^12^11 


(e^cosa^ ,-(a,+e )cos(a,+6^„) )] 


an -5 


2S' 


53 


'12 ^"3 2 


l'^12^ 


2£ 


att ^6f-®22V 

22 44 




Now choose 


1 ^12 ^11 


'33 


2a 


12 



q Yj_2 ^22 
2a 21 


(4 .18) 
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Define 


a = q 1^2 cos 6 ^^ 


(4.19) 


With (4.19) and (4. 18) substituted, (4,17) now becomes 

• - . « - 2 
o XX ± :> 2 A ' 53 '”5' ^”11 

2 

(4.20) 


“ 2 ^ 44 (^ 6 ^ /^22 ^ ^6 ^ ^5 ^ 6 ^ 


Since the last term can make indefinite, it must be 
eliminated . Let 

■^1 = °^(X5 xg + Xg) 

Integrating this with respect to time, we get 


- a x^ Xg 

where withnaoss of generality, the constant of integration 
is taken to be zero. If we now choose a Liapunov function 
7 given by ~ have a Y-fionction which has 

a negative semidefinite time derivative. Thus the 
modified and final Liapunov function is 

^ = \l ^5 ^ 22''4 ^53 ^5 ^44 ^6 ” “ "^5 ""6 + 

*^1 

+ 9 / f2_(5) d (4.21) 

o 

with the time derivative 


7 - -2jI;li\^3"^^22^2^4 




11 


-2? 


44 ^^6 


)Va 


22 

(4.22) 
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It is now necessa-ry to ensure that V itself remains 
positive definite. This can he guaranteed if the quadratic 
form in (4.21) remains positive definite in x, which 
requires the condition that 

33 44 " “ V4 >0 (4.23) 

The coefficients and when simplified, turn out 

to he 


^55 

and 

p 

'-44 


2 ^^11 ^dl ■* ^dl^/I-\l^^dl 


dl 


)] 


_ 1 


|(Z22 + ^d2 ” ^d2^/^ ^22^^d2 - ^d2^ ^ 


With these values of and and with a as given 

hy (4.19)? the condition (4.23) will he valid for a 
practical system and thus V will remain positive definite. 


4.3 three-machine SYSTEM 
Formulation: 


Consider now a three-machine system, wherein three 
synchronous machines are interconnected through transmission 
lines. With the same assumptions as in the two-machine 
system, the equations of motion of the rotors of the 
three machines can he written as 


M.6. + d.6. 
11 11 


'mi 


3 

.1 


E.E.Y. . sin(6 -6 ) 
1 J ID ID 


(i = 1,2,3) 


(4.24) 
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The rates of changes of the voltages and due 

to the variation of flux linkages can he derived hy 
extending the procedure followed in the two-machine 
system and we will get; 


E . 5 

^i " T^ "■ ^ii ®i ^ (4.25) 

ci 

where 


a , 


11 


1 

oi 


[1 


+ 


di 


di 


X 


] and a. . = 


(x 


xl 


11 


10 


di di 

XT^TT" 

10 01 


) 


( 4 . 26 ) 


s 

As before let e^ and 6^ (i = 1,2,3) be the post fault 
steady state values of Ei and 6^ (i = 1,2,3) respectively. 
Define 



^(3+i) ~ ''i 

^(6+1) = \ - n < i ) 


(4.27) 


With (4.27) substituted into (4.25), we get, on following 
the same line of simplification as in the two>r-machine 
system, 


'7 


*aj^j^Xr^ — a-j^ 2 C 0 2 oos 


a^^[e^ cos 


62^2 ""(^3 ©2^ cos 

^13 " ^^9 ^ ® 3 ^ ^15^ 
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^8 ~ "'^'22 ^8 ~ ^"21^®! ^21 “ ^^7 ^ ®1^ ^21^ 

~^23^®3 ^23 *” ^^9 ^23 ^ 

x^ = Xg ■■ ®'53_[^®2. ^31 ^ ■*" ®i^ ^ 31 ^ 

— a^2 L®2 *^52 *" ^^8 "*" ®2^ ^32^ (4 •28) 


No¥ define 

01 = = x^-x^ j a^ = x^-x^ ; 

= x^ ; 03 = xg ? 0g = x^ 


(4.29) 


Also define the nonlinearities 

f^(a) = (o^+e^) (03+62)1^2 sin(0^+6®2)-e^e2Y^2 *^12 

^* 2 ^— ^ ~ ^ ^ ^ ^{'^^3 ^^13 (^ 2'^^13 ^""^1^3^13 ^13 

^^(2) ~ ^ ^ 5 '^® 2 ^ ^ ®^ 6 '^® 3 ^^23 ^ ^^^"^^ 23 ^"*^2® 3^ 2 3 '^23 

f^(0) = a^2t®2®°®^12 ■" (^ 5 +® 2 ^'^°®^‘^ 1 '''^ 12 ^ ^ 

+ e 3 COs 6 ® 3 -( 0 ^+ 03 ) 003 ( 02 + 6 ^ 3 ) ] 

^5(2.) = ^"21^ ®1®®®^12 ” ^‘^4'*'®!^®°® 

+ a2^^ ©300 s6 2^ * (^^+®^)0 Os( 0^+6 23 ) ] 
fg(0) = a3^[e^cos6®3 -(0^+e^)cos(02+6|3)] ■ 

+a32 [©2 ^23'"^®^5 "*" ®2^°°®^^3’''^23^ ^ (4-30) 
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Making use of (4«27), (4«29) and (4«30)> tiie system 
differential equations (4.24) and (4.28) can be put 
in the form: 

^1 = ^4 

^2 ^ ^5 




^4 = 

_X X - 
1 ""4 

h/% ■ 

- 

« 

Xr- = 

5 

“^2 ^5 ^ 

fl/Mj ■ 

fj/Mj 

Xg = 

-X_ Xr + 
3 6 

fj/M, + fj/Mj 

X^ = 

~ ^-ll ^7 

- ^4 


Xq = 

^■22 ^8 

-^5 


Xg = 

-a^3 Xg - 

- fg 



(Note that fj^(2) is written here as f^^ (i=:l,2, . . . ,6) 
for brevity). 

Define new state variables 
^1 = ^4 

^2 = ®2 ^2 ^2 ^5 

= D_ X, + M_ Xt- (4.52) 

3 3 3,55 

"With these new variables replacing x^» ^2 

x^t the system can now be put in the standard form 

(4.12) with the vectors and matrices given by 
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"4 I 

i (£) 1 

1 1 : 

{ ; 

i o * 

1 1 

1 ^ 

I 

ri i 

X- 1 

5 1 

! f2(£) i 

j 

i 

1 ^2 ! 

! ; 

r =1 5 1 

1 ; 2 i 

^6 

i 

f f(5) =* I 

1 *^3 i 

5 5=1 '’1 

1^3 

x? 

1 ! ^4^^) 

; 1 4- 1 

i ■ - 


00 

f5(2) 

hs! 

t i 


Xg 


1 



-diag(^T , 


a. 




1 ' 2 ’ 3 ' 11 22 ' 33 


r i/Mi 

■-1/^2 

i 0 

! 

! 0 
0 

! 0 

r'l 

I -1 
! 0 


l/M^ 0 

0 l/M^ 

-1/M^ -1/M^ 

0 0 

0 0 

0 0 


1 

0 

-1 


0 

1 

-1 


0 0 
0 0 
0 0 


1 

0 

0 

0 

0 

0 


0 

1 

0 

0 

0 

0 


0 

0 

0 

0 

0 

1 

0 

c 

0 
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-M /D 


0 

0 

0 

0 


0 

M^/D 


0 

0 

0 

-M2/D2 

M3/D 


0 

0 

0 

0 

0 

i 

0 

0 

I 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

< — 1 

H 

t . . 

“1/D2 


"l 

0 1 

! 



l/Bl 

0 


-1/1)3 




1/^2 


-I/D3 

i 


0 

0 


0 



0 

0 


0 ! 



I 

0 

0 


0 





Liapunov Function for the ^-Machine System: 

Tentatively let us assume a liapunov function 

Vq as 

^0 = +’-22^5 * * S5^l + *- 66^9 

®i 

5 ^ 

^ I % / 

i=l ^ 0 "" "■ 

Bventhough we have a term in (4.34) which is the sum 
of integrals of nonlinearities which are of 'multi- 
argument’ type, the "behaviour of "V"^ a,round the origin 
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will "be similar to that of the single machine case. 
Thus it will be positive definite in a finite 
neighbourhood of the origin. Differentiating (4.34) 
along the system equations in the standard form and 
assuming 

hi = I h ' "22 “ ^" 2 / 2 ' S 3 = 

we hsi-ve 


^0 ~ "^\l’'l^4 " ^‘■33^3^6 


-2x 


- 44 


C “a,, -1 Xi-.~ 


‘^1^12^11/- 


7 ^‘11 7 2!?^ 


■(e2Cos6i2- 


ai 


44 


(ac-+e^) 003(0^+6^2) ) 


s N ^ ^^13^11 


( e^ cos6^^ 


5 '"2' - ^ -44 

(o.+e )cos(cr2+6® ))] + (qT^2C°®^12 ^^7 + 




a 


22 


'^^12®'22 


qY 


23 22 


2?. 


55 


(eiCos6®2-(G4+ei)cos(o^+6^2)) “ zz^c: ''-3'—-23 


(e„cos6) 


■55 


(Og+e^)cos (0^+6 2^ ) ) ] + (qY^2°‘^®^12)^7^8’^^'^'^23'^*^'^^23 )^9' Y 


'3 23 


I 


66r 


)) 


2^9 ^^”^33^ ^2l^^(ei'^os6®^-(o4+e^)cos(02+65;^) 

^,, 23' 3.2 ,(0^cos6^^-(q^+e2)cos ( ^ ^ ^13 ^^7^9 

^'■66 


+ (qY23COs6®^)xQXg 


(4.36) 



Now choose and I such that 




2 9.. , 


a 




12 


2 *-44 


a. 


13 ' 


^^12^.22 


2i 


a 


55 


21 ^ 


^'^23^22 

2 


a.. ; !!i2!u 


23 


a.- . ^h3^33 


■66 


Also define 


31 " 2£ 


=: 


66 


32 

(4.37) 


“l = ^ ^12 ^12 

cc^ = q cos 6®^ 

“3 = ^ ""23 ^23 


(4.38) 


With ( 4 . 37 ) and (4.38), the expression for reduces 


to 


Vo = -2P-,,X,xf - 2£..X,x^ - 2?^,a,x^ - 


'11 ^ 22 2^5 

!^r; - "^66- 


2 % . . 2 

(x„) - 


'33 3^6 


"22 


^^8^ arr^^4^ ^ ^ 7 .^ 8 ^ ■'' 

33 


CC^iXr^X^ + X^Xrj) + a^(XgXg + ^ 8 ^ 9 ) 


(4.39) 


By the same reasoning as in the two-machine problem 
seen ea,rlier in this chapter, we will ensure the 
negative semidefinite property of the derivative of 

the final Liapunov function chosen, by eliminating from 

« 

( 4 . 39 ) those terms that may make V^ indefinite. Thus 


define 
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Mow subtracting this v from we get the modified 
and final Liapunov function as 


V = + ^22^1 + h3^6 + * Ss^s + "66^ 

3 

“l^7^8 ^ *9. I / f^• (£)da. 

i=l 6 ^ ^ 


(4.41) 


with its time derivative, 


2^. 


v = — 2 o X X ^ A X ^ ? X A X ^ — 44 / „ \ 

^ ^11 l '^4 ^ 22 2^5 Q [^> 7 ) 


33 3 ^ 


2^1 


S5 


^•22 ® 


• 2 2 

(^«) -r^(ic) 

aj5 9 


( 4 . 42 ) 


It can be shown on appropriate substitution and 
simplification, that 

xi, ) 


£ 


44 


1 i.hl_h^dl. 


2 - ’'dh 


dl 


,, = a ilxlfdipiMsl 

'55 2 T. . fx - _ - x’ 


^22^^12 


, ^ a ^^3? ^d3 

'66 2 


“■ X j -y ) 


M3' 


S3^^d5 ^ ^d3^ 


d2- 

(4.43) 


Clearly V is negative semidef inite . Without loss of 
generality q can be chosen as unity. The positive 
definiteness of V can be guaranteed if the quadratic 
form in (4.41) is positive definite which can be 
ensured if the following three conditions derived by 
applying Sylvester's criterion, are satisfied; 
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(i) 

1/C^ - cos"6®2 > 0 


(ii) 

1/C^ - cos^d®^ > 0 


(iii) 

1-(C^cos^6®2+C2COs^6®^+C^cos^ 6®^. + 
C^cosd^^cosd^^cosd®^') > 0 

(4.44) 


where 


= 


~ ~ ^11 ^22 
^12 ^\l^^dl~^dl^ ^^22'^^d2“^d2^ 


C^ = 


(^dl " ^ 


dl' 


‘'dd ~ ^d^^ ^1 


^13^^1l'^^dl“^dl^ ^^33'‘’^d3~^d3' 


C 


= ^^d2 ^dP-^ ^^d3^^d^^ ^^22 ^33 

^ ^ 23 ^^ 22 '^^d 2 "^d 2 ^ ^S3'^^d3"^d3^ 


'4 


^^^dl~^dl^ ^^d2 ~ ^d2^^^d3 ~ ^d3^ ^11 "^22 ^dd 


Zn o^T 'T-^r^-zc ( Z-I 1 +Xj -, ~Xj, ) (Z_„+X . 


) ( Z™~+X , „— X T _ ) 


"12-^13 23' 11 dl "dl^ '^"22'^d2 ^d2'^ 33 d3 d3' 

(4.45) 


For a practical system, the conditions (4.44) are 
alwa.ys satisfied, hence ensuring the overall positive 
definiteness of V in the region of interest. 


Illustrative Example: 

Consider the system shown in Figure 4.1. 
Assume the values of the various reactances as 
follows (all in p.u.): 



109 



FIG. 4-1 A THRJiiE-rMACHIFE SYSTEM. 



no 


X 


'dl 

1. 

.1 ; 

^d2 = 1-15 

* ^d3 

'dl "" 

0. 

.23 ; 

x'g = 0.35 

o T 

' ^d3 


The transmissiGn line reactances 'between the ith and 
jth 'buses denoted by z . . ares 


==612 = °-2 ; ^ el 3 = ’• ^ e 23 = °.-5 

Assume a fault to occur along the transmission line 1-2 
which is cleared by opening the line 1-2,. For such a 
system, the rcractances involved for computing the 
coefficients in the Liapunov function can be verified 
to be 


^12 

3. 

62 ? 


1.065 

f 

^23 “ 

-P- 

CO 

TT 

^11 “ 

0. 

822^ 

^22 = 

1.025 

D 

f 

X,- = 0,. 
33 

.61 

Th^..sc:; values 

1 result in 

(with q = 

1). 


44 

1. 

19 ; 

'■55^ 

= 1,11- 

and ilgg = 1.35 


0. 

27700 

36^2 : 

cc^ - 

0. 

94cos6®^ i 

r 

II 

0. 

697co 







Thus the Liapunov function for this system is given by 
Y(x,a) =-^x| + + ^Xg t 1.19x^ + l.llXg + 


l*35x^ _ 0.277x^xgcos6j2 0. 94x^x:gCOs6^^ - 

0.697xoXqCos6® + I / 

° ^ i=l o 
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II can Toe verified that the quadra.tic part is 
always positive. This can be shown by directly applying 
Sylvester’s criterion to the above liapunov function 
or equivalently by applying conditions (4-44) wherein 
we have 

= 0.0145 ; = 0.1375 ; = 0.081 ; 

= 0 . 0256 . 

Thus the above V-function will be positive definite 
in a finite neighbourhood of the origin. Beyond this 
region, V may become negative due to the contribution 
of the integral terms. The region of. interest, however 
will lie within this neighbourhood. 

Several values of the external reactances x^. 
spread over a wide range, were tried and it was 
observed that the conditions (4.44) were always 
satisfied. Thus the V-function developed here, 
seems to be a perfectly valid liapunov f ■unction, I'/hich 
can be used for the study of all aspects of stability. 

4.4 EXTENSION TO k-MAGHINE SYSTEM 
Formulation of the Problem: 

On close observation of the two-machine and 
threo-machine system formulations, an extension to 
a k— msichinc system can be made by induction. 
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The differential equations governing the rotor 
motions and the variation ef the voltages proportional 
to the field flux linkages for the k-machine system 
are given respectively by 


M.6 +D.6 = ~ I B.E.Y. . sin(6,-6,) 

1 1 1 r mi 1 D 3 . j '1 2 

( i=l , 2 y . • ♦ ) k) 


( 4 . 46 ) 


k E 


T' 


01 


{ i“l ^ 2 ; • • • ^ k) 

A 

Define now a ^^“dimensional vector X as 


d¥i- - 

(4.47) 


Z = 


X. 


1 1 


^2 ! 


(4.48) 


■ -3 ! 


where the k-vectors X-j^, X^ 3 ^d X^ are defined by 


! ^ i 

x^ i = 6-6®; X^ 

X = , 2 1 2 

~-L 1 . i 


i 

! X- 


k 


k+1 ; 

^k+2 ! = ii ' 

i 

! 

' • ■ ! 

i ^2k J 


I X 


X. 


X 


2 k+1 
2k+2 


= E - e 


! X 


5k 


(4.49) 
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Here ^ and E are column vectors with components 

(^ii, o^j •••> 6 ^) and (E^, E^y E^) respectively. 

s 

6 and e are the post fault staady state values of 
i and B respectively. With the aid of ( 4 . 49 ) » 
the equation (4,47) can be rewritten just as in the 
two and three-machine systems, in the following form 

k 

^2k+i ^ "^-ii ^2k+i ' 1 . 

J J- 

(i*l, 2 ,...,k) (4.50) 

v.'hero 


=^11 = fhi ^di - 


13 


'01 11 ' 


(4.51) 


21 

Ijofino a kxm matrix K as 


i. = [Ki 


bc-h 


where 


^'( ( j-l)x(k-j) ) 


^3 " 


1 


(lx(k-j ) ) 


(j=:l,2, ...,k-l) (4.52) 


■^( (k-j )x(k-j) ) 


and 


m = k(k-l )/2 
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The subscripts of the form (ixj) attached to the 
various submatrices here and elsewhere denote the 
dimensions of the respective submatrices (i and j 
indicating the number of rows and columns respectively) 


In (4-52) o.bove ) Is a row matrix having all 

its elements equal to unity and 
unit matrix of order (k-j). 


■( (k-j)x{k-j ) ) 


IS a 


Now define a (m+k) -vector 

r y '! 

-1 


a = j 


I 

“2 


(4.. 53) 


whore the m-vector and the k-vector z ^ a^re defined 
by 


'1 I 


T 

IX. , 


1 ' -2 


! o ! 
m : 


m+1 ! 


m+2 i 


m+k 


= X. 


(4.54) 


Also define a {m+k)-vector f (o) , respresenting the 
nonlincarities as 


1.(2) 


%(a) i 

02(2) I 


where 

^^(a) = col(f^(£), 


^^(2) = col(f^^^(2)> 



(4.55) 


(4.56) 



115 


Tile components of ^^{o) are defined by 


fi(2) - 


(1=1,2, ., . ,in) 


(4.57) 


where j and p are the indices of the components of 
on which is dependent. 

The k-vector ^^(o) has components defined by 




37^1 


i^l ,2,...jk) 


(4.58) 


where 

°ij = “p ^ 

- -a if i > 3 
P 


(4.59) 


where o is that component of the vector which is 

dcT'endcnt on components of X^ having indices i and j 

(i.e. o =x. -X.). Using the various relationships 
p 1 3 

given above, the system equations (4.46) and (4*50) can 
be c^vSt in the follovj'ing form 


-1 ^ ^2 

= -D Xg 


M ^ K ^^(o) 


” ~-^o ~5 ^2^—^ 


(4.60) 
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where the kxk diagonal matrices M, D and are given "by 

M = diag(M^) j d = diag{D.) ; A^ = dlag(a^^) ( 4 . 6 l) 

How introducing k new variables denoted by a k-veetor^ 
and defined by 

i = D +• M (4.62) 

the syotem can be represented by the standard form 


X = F I - G f (o) 
£ = -N f (a) 

0 = h" ' X + P ^ 


(4.63) 


where Z 


-2 


; P is a 2kx2k matrix; G is a 
(2kx(ffi+k)) matrix; K is a {kx(m+k)) matrix; H is a 
(2kx(m+k)) rntrix and P is a ((m+k)xk) matrix (all real 
and constant matrices) and are given by 



°(kxlc) 


m"^k 

(kxk ) 

p = 


; G = 




1 ^^(kxk) 

1 

^ (kxm’) 

t 

^ (kxk ) 


N = [K ^(Jj;x3£)^ ^ ® 


T -1 

-K MD 0 


(mxk) 


0 


(kxk) ^(kxk) 


T -1 
KB 


0 


(kxk) 


(4.64) 
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liapuncv Ftinction for the k-Machine System; 

Usliifc the method of integration by parts as in the 
two-machine and three -machine systems., we will get a 
iinrunov fimction of the form 



where 3 and p are the column indices of the first and 
second nonzero elements respectively of the ith row of E . 
While V is seen clearly to be negative semidefinite, the 
positive- d'., finiteness of V can be guaranteed by imposing 
conditi.oiis similar to the ones given in Section which 

in essence requires that the quadratic part in the 
V-function (4..65) should be positive definite.. These 
conditions can be obtained by applying Sylvester's 
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criterion as before. It is reasonable to expect that 
in real? Stic systems, these conditions ¥ill be satisfied. 

4.5 Oe-.'CLLblOI] 

Ine theme of this chapter has been the incorpora- 
tion of flux decay effect in a multimachine system, which 
h.'ic hith-fto boon investigated with the assumption of 
ccnolr.nl li'.'.lc* flux linkages. The immediate effect of 
the addition of this factor is to render the problem more 
corr.i 1 j c'' tc fl . The methods that have been used successfully 
to derive Liapunov functions in the earlier chapters, 
could i!Ct be extended to this problem. A version of 
'int gr"tion by parts' method hr.s been found suitable 
to gerf.rato Liapunov functions for such systems. A 
two-machinc system and a three-machine system have been 
f Oi’L.ulat i.,d and appropriate Liapunov functions systematically 
dcriviui. TheS'-. results have been used to formulate the 
problcri of k-machine system and to arrive at a Y-function. 

It is interesting to note that, if the flux decay effect 
is negl; cted , the V-function derived here resembles the 

rj Q 2 ^ 

on(. rf’;y typo functions derived in the literature ’ ’ 

The voltage regulator and governor dynamics can 
in Lpl'. , be incorporated into the model examined 

in this ohept r, on similar lines adopted for the single- 
rac.chir.u system of Chapter III. It is to be hoped that the 
results of this cha.pter will form the basis for detailed 
numerical exp^.r imentation. 



CHAPTM T 


lOWiiR SYSTEMS l-ITH TRAISEER CONT'UCTMCES 

5.1 i:.T?Cj>l'CnON 

i*/.is chapter is concerned ■with, seeking Liapunov 
luuctiorio vhich can be systeiiiatically constructed for 
. c’ er tores including the effect of transfer 

c;r. oeo. Although most transient stability studies 

ill tiie cruirt have been conducted with the transfer 
con/: act- Is c '..'S neglected, in real systems these are far 
I'l-sim nor.ligibl'., in particular due to the constant 
’rr;- d-'inc’. loads that might be present at the buses. 
Ih'.u-i.i’u.; 0 rn .accurate analysis may require the 
ir. cl us lot, of the same. Inclusion of transfer conduc- 
t'.,TiCt.£i 'ii, vS not precont any difficulty in the 
a;! r '! I ;'t, : on inc.thod of stability analysis. However, 
vh'. !i til',, str bility investigation is made through 
Liciui.ov ."'Cthod, this additional detail introduces 
oi.r..-. i : ■•'■fj ciiitics. Although the number of nonlineari- 
ti..f; - .',ch of v/hich satisfies individually the 

'ci.ciuii' conditions' - becomes twice that for the 
triinsfcr conductances neglected, one 
’.vill i.np. ct that a matrix version of Kalman's 
const. auction procedure must be able to lead to a 
li-tp. t;c-v function. Indeed it does succeed in 
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yi.ldiiie a Liapunov function for a two-machine system, 
.out as v.ill be shown later in the chapter, the problems 
cl Cj .. -oiiiE with three or more machines are not amenable 
to tj. ; pi'ocoduru as the matrix Popov frequency condition 
G.' nnct bo srti&fied. Phis problem seems to have evaded 
a t ! 'Ctor.y solution in the earlier works a.lso, as 
d i on. . ■ * a t r i e f 1 y b ul o w . 

Aylott^ considered the inclusion of transfer 
c-o.fiuc '.•••■.r GOG in his energy integral formulation, 
i-w.-o V c i as pointc.d out by billiams these energy 
i r. ' 0 ..i'* -1 ■ •'■'i i,. not strictly Liapunov functions, 

^'li'l 'vu tiicsc functions can be used for stability 
studi'S,:. i'lOreovor, Aylett neglects the transfer 
com, uctr nccs in certain terms when they present 
d i !' ri cult j es in integration while retaining the sam.e 
in Cert'-, in other terms where they do not pose any 

p O 

cCi’.spiitet ii nal difficulty. Ribbens-Pavella recasts 
the tmc.blcifi a set of first order differential 
-. 00 --I.,; era: in the k(k-l) dimensional phase space 
(v;. .rG k is the number of machines) in order to 
nci.i'.vf. v.'hrt ho terms 'symmetry'. He then proposes 
a function making assumptions analogous to 

ii.-.nt, c,f Aylt.,tt. The procedures followed by both 
thi autiiors rxe applicable only for uniformly damped 
0 '/:■ . (In fact they neglect damping). Further there 

is Imclr. of consistency since the transfer admittances 



121 


are asrsuiacd as purely susceptive in some terms while in 

seine other terms involving the same pair of indices no 

such assumption is made. In the process, some terms 

with indefinite sign are ignored in V which is 

ultimately shown to he identically zero^'^. El~Ahiad 
7 

and Fast-apaan have included transfer conductances in 

their Liapunov function but have rightly restricted 

its applicability to Lagrange stability which implies 

boundedness of trajectories. Ltiders^^ proposes a 

Liopuinov function for such systems, again with some 

approximations (neglecting the nonintegrable part of V) . 

Im. ro suiting V- function however, is practically the 

sMn<.' ".3 the one he proposes for the system without 

tr^naf r I’onductances (except for a constant multiplier 

nging in value between 0.96 to 1.0, associated with 

coKn., l .cmc in V). Doubts have been raised whether the 

/-funct.ioi. d^^rived by liidors is not a ’first integral’ 

64 

jnnh'ui of : Li^'.punov function 

it is therefore evident th'^'t the inclusion of 

ti-u..','.' ccnduc bant os is not straightforward. A two- 

■nAn;'' sy..-twr cvf.n with damping included, can however 

25 

bi- h’-ia'led. Lai cand Mui-thy ha.ve derived Liapunov 

I'l..; ai i-!u ruch system using a generaLization of 

35 

ior'ov’;' cfit-i'ion due to Moore and Anderson . This 
’ .'iO I, ; r I! h Ir.vcKtigpt^d in this work through 
halm- a'.; cc..j t .-uction procedure. The problem of 
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three-machine system will he formiilated to exhibit the 
'lack of symmetry' in the transfer function matrix, 
xvhich precisely is the reason why the Popov frequency 
criterion is not satisfied, resulting in the failure 
of this particular line of attack. 


5 . 2 I't/’iD-MACHIhE SYSTEM 

Formulation of the Problem: 

With the usual assumptions including constant 
field flux linkages and constant input and neglecting 
tronsi,,.nt saliency the two-machine system with transfer 
conductances can be characterized by the following 
dynamic equations 


+ D, 


- hi 


+ Dj 

®2 = ^m2 ■ 

■ ^e2 

(5.1) 


The- electrical powers P^^ and P^^ given by 


Be Et E« 

1 . 12 . ( c 

■ =,, — Sin a-,-, +-U 5 sin(o 

11 


12 


' = ^ sin n ^ 




“12^ 


“l2> 


( 5 ^ 2 ) 


W!-;.:., r<- now Z,. .. and Z. . are the driving point and 

.X. 1. ,3 

i:,:pud‘ rices rc.spectively y corresponding to 
tLi, tv* perf omiB-nce of the synchronous machines 

(if !.i: L»nit iir.T C'd*inces) • and oc^j are the compliments 



of -n'l respectively; and are the 

i.v,-r:c. maos of Z-^ and Z-^ respectively. Farther 
7 - 7 md 6. . = 6.-6., Substituting (5.2)into 

^ij ■" ji ' ^ ^ 

(■^.i; and simplifying, we have 


T' E B 

h 6 - i- [4-2 sin (6,, -a,,) 

% 1 % Zi2 12 12 


!^sln(6^, - «i2a 

h2 


“2 6 - i_[ ^ since,, 
■ 2 Wj h2 


%®2 


sin(6®i 


Ui^)] 


( 5 . 3 ) 


Note that in arriving at (5.3). -ti^e term 

(p . - a sin a., ) has been replaced by the secona 

^ii ^ . nf fS 3), in terms of th< 

term in the square bracket of (5- ^ 

post fault static e<iullibrlujr. points .6^ (i=1.2 • 
Define 

cS 

X = 6i - 


c ^ 

Xp = 6, - 6, 


yi = h = h 

• * 

yg = ^2 " ^2 


(5.4) 


Define also 


Oi = ^1 - 

Op = - x^ 


(5.5) 



i'jt tf’..', ncnlinearitles be defined as 


^ -■ “12) - - c<i 2 )] 


ii t E 

= 2 ^ [sin(a2 + 6|^- - 0i^2^ slni^l 


21' *■ “12^ ^ 
(5.6) 


Using equations (5t4) to (5-6) and with 

and ^2 ~ system equations (5-3) can he 

c'’st in the form 




X 2 = y 2 

^ 1^1 " 

^2 " “ ^ 2^2 ” Mg 


(5.7) 


with o^ = x^-Xg and Og = Xg-x^. 


Define two new variables ? ^ and Cg as 


?! D^x^ + M^y^ 
?2 = ^ 2^2 ^ ^^ 2^2 


(5.8) 


With these new variables, the two-ma, chine system 
( 5 . 7 ) can be recast in the standard form as 
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1 

1 
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^2 
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IkI 

0 

i 

1 1 
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1 
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The tr-insfcr function matrix of the linear part of the 
above system can be derived to be 


I 

; s ( s+ ^ ^ ) 

I 

¥(s) » I 

I _ 

I s'(s+'^) 

L. 


1/Mg 

s(s+X2) 

1/Mg 

s(3+^^^) j 


( 5 . 10 ) 


Liapunov Function for the Two-Machine System; 

Now for the system (5.9) a Liapunov function is 
to be derived. This requires at the outset a choice of 
matrices R and Q such that the following matrix Popov 
frequency condition^ ’ is satisfied; 

Z(3 w) + Z^(-jw) > 0 (5.11) 


where 

Z(s) = [2RPN + Qs]¥(s) 


(5.12) 


In (5.12), the matrix R must be symmetric and so chosen 
that the product RPN must be a diagonal matrix. In the 



rr'.'cl :r. undi.r study, the matrices P and N are such 
tnc above condition can be satisfied only by 


r.akin;; the product RPK a zero matrix, This can be 
achieved by choosing K as a matrix with all elements 
equal. Tv;o possibilities exist. Either a finite value 
C':n be chosen for the elements of B, or R can be chosen 
as a zero rriatriz. It is verified that either choice 
results in tho same Liapunov function because of the 
S] ci'-l nature of H and o, ¥e may therefore, as well 
chosG R as zero matrix. Thus choose 




Then the matrix Z(s) is seen to be 


( 5 . 13 ) 


Z(s) 


S 4-X^ 
s+Xjj^ s+Xg 


S 4- ^ /-v 


^ 5/^2 


( 5 . 14 ) 


The matrix Popov frequency condition then requires 

ZX^qi/Ml 




<ll X 1 ^2^ 




M 




l-i- ] 

M. 


M 


1 


X + m * 


• • ^1^ 1 ^2 ^2 






1 


> 0 


''M2 ^ ^M2 

i 

jw^4-X^X2 “ ^ ^2^ 


M. 




+ 0,= 


(■5.15) 
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.r.i;- i'./r.i.rc£j tht; restrictions 

‘ 2 ^ ^ “ ^ 2^2 (5.16) 

Tri.. ...lU' s;ions (5.16) imply 

1-^2 (5.17) 

Thu;.; it is ncces&'^.ry to assume 'uniform damping’ 
fcr ti'i."; ryst^^m if v'o are to succeed in arriving at a 
li-''*pur.ov function. Furthermore, one of the elements 
in uhc diagonal of Q can he chosen arbitrarily , while 
thw othL-r clement is fixed by the relation (5.16). 

V/ith R and Q chosen and having established that 
the matrix I-opov frequency condition is satisfied for 
the two-machine system with uniform damping, we proceed 
to construct a Liapunov function for the problem. The 
procedure to derive the L matrix is along the linos 
outlined in Chapter III through equations (3*25) ** 
(3*28)'. However, because of the special nature of the 
matrices involved and the assimiption of uniform damping, 
there is no need to do explicit spectral fs,ct orization ■ 
involved in the procedure. The general technique of 
handling this problem in a .k-maciiine system has been 
outlined by Pai •. A Liapunov function for the problem 
on hand is now derived as follows: 

(i) A 2x2 matrix T(s) is to be found such that 
i[Z(s) + Z^(-s)] = T^(-s) T(s) 
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With P-nd 12/^1 == (5.16), it 

c " i". b * V r i f i u d t t 

i[2(s) + Z^C-s)] — T (5.18) 

(s+X)(-s+x) 


j *11/^1 ■''ll /“2 

' I (5.19) 


Thus 

T(s) = T (5.20) 

S +X 

Although in general T is to he determined by 
factorizing (5.19), iu this case, it is possible to 
avoid the explicit factorization as can be seen in 
thu sequel. 

(ii) The matrix ^ is next determined from 

+ PN) + H+ P^)Q^ = r 

With the appropriate matrices in the above problem, 
it can be seen that 

r = 0 ' (5.21) 

(iii) The ma.trix U is solved in general from 

T(s) _ r * -U^(8l - G 

which on substitution of the various matrices yields 
/X T = G ( 5 . 22 ) 
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, m 

i. jr, ' 1 -explicitly solving for U , we proceed 

fiij',o‘l;; to oclvo the rsotrix liopunov equation. 


Ov; The Liapunov natrix equation is given by 

P^l + I? = -U (5.23) 

T 

i r*. ti plying and post-multiplying (5.25) by & and 
1 r-et.poctivtly , we have 

I + iP)G = -G^UTJ^G (5.24) 

yith P = - A I and with (5.22), equation (5.24) 
ri-ducoP3 to 

m ni V 

2G^LG = T T (5.25) 


Assuine the symetric L matrix to be of the general 


forri 


L 


r 




12 

22 I 


( 5 . 26 ) 


Substituting this in (5.25) and solving for the 
elements of the i matrix we get 

~ *^1^1^^ ^ 12 ~ ""'^l^l-^^ ' 

^22 = ^ 2 ^ 2 /^ = U«l/2 ' '' 

The Liapunov function is therefore given by 

2 

T = 1 y + I oi /"■ fi(Oi>aOi 

1=1 0 

which on substitution of the above ma.tr ix 1 and 
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jr.tnt 


yields 


t f ^ 


•'iCyi-yp)^ + 
1 E^E, 

. / 

^ • "12 


*■ 1 * 1 '•‘^1 •>’ 2 ' 
c 

X Ju 2^ 

p}]dcr^ + 


, sin(a^+6^^-a^^)-sin(S^^--aj^^) ]da^ 

(5.28) 


c "12 


w i th 


? 


■<iiMi ^(yi-y2> 


(5.29) 


Tiu; ¥-functicn derived above car be used to study all 
aspects of transient stability of the t¥0-Eiachine systen. 
This Liapunov function differs froia those in the 
literature for two-mo. chine systems in that the quadratic 
part does not represent the kinetic energy. 

5 . 3 FORMU’LATION OF 3-MACHINE SYSTEM 


The equations of notion of the three-machines 
car be written as 

Vl+Vi = hi"^el (5-30) 


The electrical power outputs of the machines are given 


by 


ei 



sina 


11 



E . E . 

sin(6^j-a.p 


(i=l,2,3) 


(5.31) 
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cti? h'-re that Z, . = Z.. and ^ • Substituting 

•* J JX X J JX 

(b.31i into (5 -30), we see on post fault steady state 



3 ii.E. 

Gin a,,. = y -X—l sin ( 6 ® . 

'‘ii A 

dVi 

(i = 1,2,3) 



(5.32) 


Uoing- tije above equations, the system (5.30) can be 
rovritten as 


6 i 


■ Hh - r J , 

X iJ “"X. X J 


jVi 


(i=l,2,3) 


(5.33) 


Define now, 

-s- — /) —6® • T — 6 —6^ i ^'5 ^ ~ 

^1 ” "^1 °1 ' ^2 " 2 ^ ^ " 

yi = ^1 ' ^2 = ^2 ^ ^ 3 ^ ^3 


(5.34) 


Also define 


Oi = X 3 _>-X 2 ; ''l '''3 ’ ^3 “ '" 2"''3 


S 2 e: X — X-, « 0 *-- ~ 

4 2 15 51 o 94: 


(5.35) 


Let 

^ ^2 ^ ®1®3^^13 ' S ^ ®2^3’^^23 

(5.36) 
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T i* : 

,i,<( ‘ ’ 


Til., ti;o r.onlinearities as 


f ^ ( 0 ^ [ h in ( c^2‘^^12“^12 ^ ® ^ ^ ^ 

fp(c^) = Kp[sin(02+^^^-a2.5^ ~ sin(6®^-a^ ^)] 


13 13 ' 


ly.(o^) ~ ix^Csin ( sin(6^.;^—a^^) ] 


3 23 23' 

.s 


23 23' 

s 


^ 4 (^ 4 ) = sin ( 0 ^+ 6 2 ^- 02 ^) - sin(62^-'a2i)3 
fc5(0^) = K2[sin(0^+6®^-a^^) - sin(6®j^~a^^ ) ], 

fg(og.) = Z^[sin(0g+6®2“*^32^ “ sin(6®2'«^2^^ 

Using equations (5.34) to (5.37)» the system (5*33) 
can he written as 


y. 


"3 = 


■^1^1 


■^2^2 


■-K^-x 
3 3 


■ ^3^^3^ 

■'M“ ^5 ^^5^ 






SI ^6 ^'^6^ 
5 


( 5 . 38 ) 


Define new variables 


5^ = D^x^ + M^y. (i=l,2,3) ( 5 . 39 ) 

With these new variables defined in (5.59) > the three- 
machine system can be recast in the standard form (1.8) 
with the various vectors given by 



~ c ol { 1 
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( 5 . 40 ) 



( 5 . 41 ) 



154 


■„ ith ^.h>, trices defined as above, the transfer 
fin.ction rrtrix of the linear part of the above system 

t;i,ver. by 

Vl(s) = H^(sl - G + i P F 

s 

is found to be 



Wi(s) w^(s) 

-w^ ( s ) 

-W2(s) 

0 

0 

■i 


w^Cs) w^(s) 

0 

0 

-w^ ( s ) 

-w^(s ) 


0 0 

vr^(s) 

w^Cs) 

-w,{s) 

-W,(s) 1 

W(c)= 

-w^(s) -W^(s) 

vr^(s) 

W2(s) 

0 

1 

0 i 

i 


-•W^(s) -w^(s) 

0 

0 

w^(s) 

w^(s) 1 

j 

0 0 

-w^Cs) 

-w^ ( s ) 

w^(s) 

w^(s) 1 






(5.42) 

where 






W^(s) = 


v^(s) = 

1/^2 



= ‘ > 

s(s+?u ) 

s (s+X^) 

T 


w^(s) = 

1/M, 

s (s+X,) 




(5.45) 


Inspection of ¥(s) reveals that the transfer f miction 
matrix for the three-machine system does not have 
symmetric nonzero pattern, hith this ¥(s)j the matrix 
Popov frequency criterion is not satisfied even under 
the assumption of uniform damping. Hence the method 
fails in arriving at a Liapunov function. However, 



155 


r ' cauisi- ci thfi sufficiency nature of the Liapunov 
thocrt,"., *'.o conclusions can he made about the stability 
nr ino'G ihilify of the system. A similar outcome is 
oh 'u rv ; I'l systems in which more than three machines 
at'', iTcront. It can be easily seen that if transfer 
c one u c hone cs are neglected, then 

Ii{Oi) = (i = l,2,...,m) 

ar (1 VO need to define only m nonlinearities. The 
r-robleiii then reduces to the standard case which has 
bi.cn ’.f:'': octivcly solved in the literature. 

5.4 COhCLDSIOh 

In this chapter, the power system stability 
has been investigated taking into, considers.! ion the 
transfer conductances. In goneral, it is observed thet 
the number of nonlinoarities for such a system is 
twice that for the corresponding system with transfer 
conductances neglected. A two— machine system hsis been 
formulated and a Liapunov function derived through 
Kalman's construction procedure, although 'uniform 
damping' hs.d to be 3 .ssumed . Such an as sumption, however 
is not unusual, as the inaccuracy involved is snmll. 

The three-machine system has been formulated on similar 
lines as the two-machine system, but a Liapunov 
function could not be constructed because of the 
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luru rf the 'Fopov frequency criterion V' test . 
;,rublcn calls for further work as no systematic 
■•.cceirctc method to construct a liapunov function 
:i,vEtuias with more machines than two, has yet 


li found. 



CHAPTER VI 


CONCIUSIOiV 

6.1 IHTKOUTCTION 

Interest in the application of liapunov method 
for the transient stability analysis has been gaining 
momentum in recent years. After an early breakthrough 

7 

in this a,rea by El-Abiad and Ragappan followed by 
rmny others, there seemed to have followed a lean 
period v?ith regard to the further advance in the 
ai'pllcation of this approach. In fact there appears 
to be some degree of pessimism on the part of the 
power industry regarding the utility of this method. 
This can be attributed to several factors. One of the 
striking drawbacks of this approach is the nature of 
the conservative results obtained. V'hile this ma,y 
be pa^rtly due to the sufficiency nature inherent in 
the Liapunov theorem itself, it is aggravated by the 
approximations made in modelling the system and the 
consequent imperfections in the liapunov functions 
developed. In this context, the synchronous machine 
is the most vital element wrhose exact modelling is 
quite complex. Factors like damping sciliency, flux 
decay, voltage regulator and governor dynamics are 
known to affect the stability property of a power 
system. It is therefore logical to expect, that a 
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Li-’j/.i../.' !v;r;ction derived for a model which includes all 
tii.G-.,, "/’ilr:, will yield more accurate results. Inclusion 

; f -'I] there factors in the development of a liapunov 
frrrti /U hov.*u.vcr, is far from easy. Addition of greater 
ilr in the- machine model results not only in a 
highf r crfl„-r system, but often in more complicated and 
gru';,t'-r nur;-bor of nonlinearities. An important step in 
f --rmul'tting these problems lies in properly defining 
the v.riables and nonlinearities. Certain models 

result in ’aultiargument ' nonlinearities. In such 
situations it is not always easy to ensure that the 
nonlinoaritics shall, individually, satisfy the 'sector 
cc.nditions' away from the origin, although at the origin 
they can be forced to vanish. One has to then rely 
more on the quadratic part of the Liapunov function 
dominating the contribution due to the integral terms 
even if the latter becomes negative in some region 
around the origin, so that the overall V-iunction still 
remains positive in a sufficiently large region around 
the origin, within which the stability region is 
contained. The sign definite property of the time- 
derivative of the V-function also sometimes calls for 
some skilful manipulations and assumptions. 
Notwithstanding these drawbacks, the Liapuinov approach 
could still play a complementary role in assessing 
transient stability by cutting down on the number of 
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•i.t-.Iic'I siuiiln.ticns that have to be carried cut by 
ti; . direct simlat ion method. This complementary 
r 1',, \t Liapunov approach has not been well- 
:r*,C’ ,*ni,sor' by the power industry. The recent studies 
of ‘iillini Hibbens-Pavella and others^^'^^ on 
some re^'l systems have amply demonstrated the utility 
and feasibility of the Liapunov approach. These 
researchers have appLied the method for the transient 
st-'bility study of fairly large systems and the 
results h‘-'ve been very satisfactory. With the 
incorporation of improved system models and the 
consequent improvement of Liapunov functions, the 
Liapunov approach might be more attractive and 
readily acceptable. The aim of this thesis has been 
therefore, to incorporate improved synchronous 
machine models while constructing Liapunov functions 
for both single and multimachine systems. In 
addition, the utility of the Liapunov approach in 
carrying out a comprehensive pa.rameter analysis is 

а. lso emphasiiZed in the thesis. 

б . 2 REVIEW 

With the above objective in view, saliency 
effect has been used in the machine model in the 
second chapter. Both single machine- inf inite bus 
system and a two-machine system have been formulated 
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r.r.'i r err irtt- lie pun ot functions derived. A detailed 

prr' U'-t . r •••.nr.’lyGis of the single machine- infinite bus 
f-.r. iv-monst rotes how the various parameters affect 
ti.e ir' nni nt stability properties; The results 

rly indicate that saliency in general improves the 
r: t-'tbility property of the system. The extent of the 
iriprovemont however, is also related in some way to 
other fv, ctors like the damping in the system, the 
iriiti-'.l load on the generator, the location of the 
fault and inertia of the machine; The increase in 
the critical clearing time is particularly pronounced 
with a lox-jer initial load and higher damping coeffi- 
cients. Similarly if the fault is away from the buses, 
soliency has greater influence on t than with a 
fault close to the buses. The parameter analysis in 
addition to demonstrating the effect of saliency on 
the stability behaviour, also illustrates the effects of 
other parameters on transient sta,bility. Such a detailed 
analysis by the direct simulation method would have 
resulted in excessive computational effort. A similar 
study can be made for the two— ma-chine system for which 
Liapunov functions have been derived for both uniform 
and nonuniforin damping cases; One of the specia.l 
features of the Liapunov functions that ha^ve been 
derived for both single— machine system as well as 
two-machine system is that the damping coefficient 
explicitly appears in the V-f unctions i It is 
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lc.„Lc-l thtrefore to expect that these liapmov 
:';:r;ction 3 are wore accurate than those which do 
r.oi cc rtain the damping coefficient explicitly. The 
ca.-,e of systems with more number of machines having 
o'ilitincy however, becomes highly complex and derivation 
of Lir punov function becomes difficult. This is 
because cf the fact that the resulting system equations 
contain in addition to the nonlinear differential 
oqur’tions, a set of nonlinear algebraic equations 
cont lir-.ing the state variables implicitly. Adequate 
ruGults ore not available in stability theory at the 
present time to handle such problems. The extension 
of the method, to a general k-machine system thus, 
seem to be impossible unless some approximations are 
made. One possible approxim,ation that can be made is 
to assume the voltage (i = l,2,...,k) as equal 

to the voltage (i = l,2,...,k) (the latter being 
the volta,ge behind the direct axis tra,iisi6nt reactance) 
V7hile including the second harmonic term due to 
saliency effect in the electrical power output of the 
machine. There is however, some inconsistency in this 
approximation; for, by assuming as equal to Bj_, 

it is implied that x^ = and in such an event, the 
second harm-onic term in the electrical power will 
also disappeor. The problem then reduces to one 
without saliency. Therefore this problem is still 
open for further investigation. 
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The study of single machine-infinite bus system 
:r e rn continued further in the third chapter. As a 
prelude tc a generalized model of a single- machine 
cyuuuirj, flux decay and voltage regulator action have 
tc'-.u included in the modelling of the synchronous 
machine-. The problem turns out to be one having - 
mul tiai-gument type nonlinea.rity . A straight applica- 
tion of Kalman’s construction procedure does not work 
for this case. The incorporation of results from the 
work of Desoor and Wu"^*^ in the Kalman's construction 
pi-ocodure, enables us to arrive at a Liapunov function. 
Using this Liapunov function 3,n illustrative example 
has been studied. The effects of flux decay and the 
voltage regulator on the critical clearing time have 
been clercrly shown. The results show that flux decay 
in general reduces the critical clearing time, or 
equivalently affects the transient stability adversely. 
On the other hand voltage regulator improves the 
stability property, by neutra.lizing the adverse effect 
of flux decay. The extent of neutralization and 
improvement of the critical clearing time however, is 
seen to vary with the main parameters K^ and T^ of 
the voltage regulator, the former being more effective 
than the latter. In fact, decreasing T^ (or 
equivalently increasing the speed of response of the 
voltage regulator), results in only a small increase 



ii. th' critic-^il clearing time. An augmented liapunov 
i ’metier, haa also been proposed, which gives slightly 
C'.ltcr results then the one referred to above. 

How i,,-vc„r , the choice of the coefficients of the 
ad.iitionul terms in the quadratic part of the 
V-function requires some amount of search,. It is not 
known v/hether an optimum pair of coefficients for one 
system will be good for another. This has to be 
established only after some- experimentation with 
different systems. The remaining part of the third 
chapter has been devoted to the generalized model of 
the single-machine system, wherein saliency and 
governor dynamics have also been included in addition 
to the flux decay and voltage regulator. A V-function 
has been derived for this model using the modified 
Kalman's construction procedure. 

With, the success achieved with single-machine 
system, the next logical step is to extend the theory 
to multimachine systems. Since the saliency effect in 
a multimachine system ma.de the. problem more difficult 
to handle, only. flux decay has been considered. The 
nature of the multiargument nonlinearities become more 
complex than in the single machine case. The procedure 
based on Popov's theorem were found to be not applicable. 
A version of integration by parts however, has been found 
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suit -tic. Using this technique, liapunov fimctions 
lirvc Gccn systematically derived for the two-^machine 
and ti.rcC-cia chine systems. The results have been 
tj-u n '■■y tended to a general k-machine system. In this 
r 'rn.ulrition, constraints on the parameters of the 
c.ystcr. in order to ensure the positive definiteness 
-f the quadratic part of the liapunov function, are 
found necessary. V/hile this may appear to be restric- 
tivi , it is found by taking certain typical values 
fir a large number of cases in three-machine systems, 
that these constraints are indeed satisfied. The 
r..sults of this chapter should form a basis to 
include other details such as voltage regulator and 
governor dynamics. In fact the addition of governor 
dynamics (assumed linear) poses no problem. Augmen- 
ting the V-function by the addition of a square term 
of the state variable defined for the governor 
dynamics as in Section Chapter III, for each 

machine, is all that is required. The coefficient 
associated with this term can be so chosen as to 
ensure the negative semidefinite property of V. In 
principle, inclusion of voltage regulator must also 
be possible. 

The fifth chapter has been devoted to the 
problem of power system vjith transfer conductances. 
Inclusion of transfer conductances in a liapunov 
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i.as been sucessfully demonstrated using 
i';j construction procedure for the two-machine 
f;. although limited to the tiniformly damped 

: y fL . hx ten si on to higher order systems has not 

1*-' n pur-uyible since for such systems the Popov 
frc; 'uoncy condition is not satisfied. This problem 
rom-iliis still very muc|i open for further work. Even 
in tlu; two machine system that has been solved, 
assumption of uniform damping is one limitation, 
althou.yh this is not too serious since the 
difference between the results for power systems 
with uniform and nonuniform, damping is small. 

6.5 AREAS FOR FURTHER WORK 

Based on this research work, there are still 
some unsolved problems in this field and are briefly 
explained below; 

1. Inclusion of saliency for a general k-ma chine 
system has to be investigated. As mentioned elsewhere, 
formulation of the general k-machine system results in 
a set of implicit nonlinear differential eq_uations 
and nonlinear algebraic equations and presently no 
results are available in sta,bility theory to 
investigate such a problem. Some judicious approxi- 
mations have to be made if existing theoretical 
results are to be applicable. 
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(2) orcact rupresentation of voltage regulator 
arc dc-sirablc. The voltage regulator action 

:Lo rcr.r; LBivG to terminal voltage. But, if terminal 
"It'-.gc is tc be brought into the picture in a detailed 
rc; ot,i,cntr<ti''n of voltage regulator (as opposed to a 
f.;,i mulo rcprcst-ntation as in this thesis), one again 
cnc.untors the problem as in (1) above, namely a set 
of nonlinear differential equations and nonlinear 
algebraic equations. This calls for some basic 
resvnxch in control theory. 

(3) Inclusion of voltage regulator action and 
governor dynamics in the formulation of multimachine 
systems discussed in Chapter IV, has to be persued. 
Analytically this problem may sound formida.ble, but 
the success achieved in the case of single-machine 
system in Chapter III encourages one to extend the 
results to multimachine systems. 

(4) Inclusion of transfer conductances for the 
general k-machine system (k > 3) can also be 
investigated further. This problem ha,s eluded several 
reseanchers in the pa,st. As shewn in this thesis, 

Popov based methods are not successful. It is perhaps 
interesting to see if integration by pa.rts method 
will yield some appropriate liapuncv function. 

(5) In the Liapunov function constructed in 

this thesis, for systems involving multinonlinearities 
(both single-machine with flux decay and multimachine 
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j ; it net been possible to include damping 
^:;v .] icitl” i,r. the liapunov functions. This is because 
f ti'.c inability to find a nonzero R in the multiplier 
(2 /u 1 + Qc) appearing in the matrix Popov ■ frequency 
cri t^ri •■-r:. On a close examination, it was found that 
this d.ifficulty has been due to nonsingular nature of 
I'h 'ill the consequent inability to select a symmetric 
it r.ati'ix such tho.t EPh is diagonal. The singularity of 
PH sO'^ris to be inherent in the power system models 
that have been used for multimachine systems. On the 
other hand, in the formulation based on the theorem 
of Moore and Inderson^^, the multiplier is (A+Bs) with 
A > 0, B > 0, (A+B) > 0, and A and B being diagonal. 

Here A and B are not directly dependent on the system 
equations and it is possible to include damping 
parameter explicitly in the V-functicn through this 
approach^^. The fact that the tvfo procedures do not 
yield identical results have to be reconciled. A 
modification of multincnlinear M-K-Y lemma appears to 
be necessary to make the multiplier independent of 
the system equations. This is a problem for further 
research in nonlinear stability theory. Such a 
solution will pave the way to automatically include 
the damping constants in the V-function. 
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APFETOIX A 


IliFUIOV STiBILITY 

CciiBid c'.T an autonomous system^"^' 

i - liz) ; 1(0) = 0 (A.i) 

I^t 0 |^ denote the spherical region || x||< R and 
Hj^ the boundary sphere of the region,. Assume that 
the syKtem (A.l) is of class in a certain region 

r""! 

O • 

P 

Liapunov defines the origin as: 

'Stable' for (A.l) whenever given any 0<e < p 

there corresponds to it a 0 r)(£)h. e such that if 

x(t) is a solution whose initial position = x(0) 

lies in S^ then xit) lies in S ever afters 
T) 

'Asymptotically stable' whenever the origin 
is stable, and furthermore for some e every solution 
x(t) as above ->■ 0 as t->- ; 

'Unstable' whenever given any 0<e < p and no 
matter what 0< b < e, there is always an x(t) as 
above reaching at some t > 0. 

I. Stability Theorem: Whenever for some S ^ there 
exists a positive definite fxmetion Y{x) whose 
derivative V along the paths of (A.l) is nega.tive 
semidefinite in S then the origin is stable. 
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7 ’.I't c tic Stability Theorem! Whenever Y is 

; cti 1 l.y nc£;''tive definite in S , the origin is- 

P 

•• syr.ptcticr.lly stable. 

ill. ?hc Barba shin-Krassovskii gomplement to the 
A. !;/r_ trtic St'^bility Theorem! In the theorem II 
'‘hove* let 

(■0 S. = S ; 

(b) Y(x) ■> «> with ||xjl -»■ “ 

Then the system is asymptotiCcally stable in the 
l.'irge or completely stable. 

IV. Theorem of LaSalle^®: let V(x) be a scalar- 
function with continuous first partia,ls satisfying 

1) Y(x) > 0 for all x ^ 0 ; V(0) =0. 

2) Y(x) <_ 0 for all x 
• 5) V(x) “ as 11 xll -5- 

If V is not identically zero along any solution other 
than the origin, then the system (A*l) is completely 
sta.ble. 

Remark; In theorem lY the region ^ = Soo. If instead 
we have p<°° , then the theorem is valid in the finite 

region Sp , which is termed as the region of attraction. 
Indeed pov^er systems do belong to this category having 
only finite stability domain. It is this version of the 
theorem on which the results of the thesis have Doen 


based . 



APrMDIX B 


ifiOtP GF THE THEOSBM IF SECTIOF 3.2.2 

ISEuise the Lie,p\mov fimction 

V( 3 £,c) = I + (£ " 2:)^ ^(2. ~ y) + V^(£) 

(B.l) 

Th .. time doriv^.tive of V along the system equations 
c bo shown to be 

V = y^(P^L + LF) y - 2y^(LG-HRpd^ - iF^HQ^)f(o) 
[iQ(H^G +P d^) + i(G%+d P^)Q^]f(o) 


-20^ RP d^f(p) (S' 2) 

Define 

F^L + IF = - e (S' 5) 

L G - H R £ d^- H = P (BM) 

iQ(H^G+P d^) + i(G^H+d /)Q^ = ^ (Bs5) 


With the above defined quantities 
Y = -[y^e y + 2 y^u f(o) + f^(o)K f(o)] 

-2 R P.d^ f(o) (B.6) 

Let the qup.ntity in the square brackets be S. Together 
With the assumption a2 s it is required that S >.0 
to make V< 0, Further it is necessary^® that k: >.0. 
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It. i;; th'.5C poscible to write 


K= r (B.7) 

v.-h,r*, r io !x (mxm) matrix. Substituting for k and 
with r >0, S can be written as 

L = !! r f (0) + y P - y^(e - U U^) y (B,8) 


v;hi..r c 

u r = p 

Du fine 

= 0 - U 

To ensure S> 0, it is now necessary and sufficient 
that Thus if there exist a positive &t,mi— 

definite matrix 9^, a symmetric positive definite 
matrix 1 and a constant matrix U such that U ^ = u 
and 

(a) L + I F > -(^i + n 

(b) I & - iP^ H - H R£ = ur (B.ll) 

then S > 0. The two equations in (B.ll) are the 
so-called Lur'e’s resolving prelimit equations for 
the multinonlinear case. If we take T = 0, b can 
be made positive semidefinite in y only by setting 
u = 0 and S. would then loecoiiie 


(B.9) 

(B.IO) 


T 

s = y © y 


(B-12) 
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j-i.-. ccri'u jpoiiJ ing Lur’e’s PGsolving emulations are 
::.=.r'ay r, particular case of (B.ll). Hence (B.ll) 
ccv..r>,i, in gencrol, the case r> 0. Prom (B.11) 
the .uiu’'i.;’c resolving limit equations are obtained 
i.,ttirig 9^ = 0. 

How consider the matrix version of Ealman- 
YacubovaLch lemirB^^’ 

Let ^ be a real (mxm) matrix^ and P, a real 
( (r.-l )7„(n-l) ) stable matrix. Let G and K be two 
( ( n-l ) xi:» ) m.atricos . If 

J(s) = (sI-P)~^ G (B.13) 

the frequency condition 

+ J(D“) + >_ 0 for all realm (B.14) 

is nocessa,ry and sufficient condition for the existence 
of a ((n-1) X (n-1)) real symetric matrix 1 >.0 and 
a ((n-l)xni) real matrix U such that 

(a) P^ L + L P = -U 

(b) L G - K = U r- (B.15) 

Compa,ring (B.15) with the limit equations correspoxiding 
to (B.ll) (i.e. with = O), the matrices P, G,D and 
r in both sets of equations are identical ano further— 

more, we equate 

m m T 

K = i P H Q + H R p d 


(B.16) 
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..it;. l quiv'' lonces , (B,13) becomes after some 

"1 i c irripulations, 

H^G + (Rpd^ + iQs)H^(sI-F)''^G (B.17) 

bmin,;, (3.7) and (B.17), the left-hand side of 

ti.. fr'i.'iuency condition (B.14) in the 1-Y lemma can 

b\, nhcswn to be 



= i[Z(j‘^) + Z^(-J‘^)] 

(since R P d^ P. seen to be symmetric). 

Thus the frequency conditions (B .14 ) o-nd (3-22) 
are equivalent. This completes the proof. 
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